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Preface 


The  objective  of  this  thesis  was  to  develop  a method  to 
identify  the  type  of  modulator-demodulator  (modem)  used  on 
various  communications  links  for  the  Automated  Tech  Control 
(ATEC) . I chose  to  model  and  characterize  seven  modems, 
choosing  as  the  characterization  the  power  spectral  density 
and  autocorrelation  function  of  the  output  signal  from  each 
of  the  seven  modems.  I also  used  the  equations  for  the 
power  spectral  density  of  five  of  the  modems  to  determine 
the  data  rate  at  which  these  modems  were  operating. 

I would  like  to  express  my  appreciation  and  acknowledge 
my  indebtedness  to  Major  Joseph  W.  Carl,  Assistant  Professor 
in  the  AFIT  Electrical  Engineering  Department  and  thesis 
advisor  for  his  advice  and  guidance  during  this  investiga- 
tion . 

And  last  but  not  least,  a debt  of  gratitude  is  due  my 
wife  and  son.  May  we  once  again  renew  our  acquaintances. 

Robert  L.  Jones 
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Abstract 

A method  for  modem  modeling  and  characterization  is 
presented.  The  properties  of  Fourier  transforms  and  the 
properties  of  correlation  and  power  spectral  density  for 
periodic  and  random  functions  are  reviewed.  The  modems 
modeled  are  characterized  by  the  theoretical  calculation  of 
the  power  spectral  density  and  autocorrelation  function  of 
the  output  signal  from  the  modem  under  ideal  conditions. 

The  input  signal  to  the  modems  characterized  is  either  a 
random  pulse  train  or  a special  periodic  function.  Also, 
power  spectral  density  curves  are  constructed  to  determine 
the  data  rate  at  which  a modem  is  operating.  Though  no 
empirical  data  exists,  the  theoretical  calculations  indi- 
cate that  modems  can  be  characterized  by  the  power  spectral 
density  and  the  autocorrelation  function.  More  research, 
however,  is  needed  for  accurately  determining  the  data 
rate  at  which  a modem  is  operating. 
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MODEM  MODELING  AND  CHARACTERIZATION 


I . Introduction 


As  more  and  more  sophisticated  communications  systems 
are  developed,  one  basic  goal  still  remains  the  same:  the 
receiver  must  produce  as  its  output  a recognizable  replica 
of  the  information  that  is  transmitted  to  it  from  a source 
through  some  type  of  medium  commonly  called  a channel.  A 
typical  system  diagram  appears  in  Figure  1. 


Figure  1.  Block  Diagram  of  a Basic  Communications 
System  (Ref  9s  2) 

Because  the  information  being  transmitted  usually 
occupies  baseband  frequencies  near  dc,  it  is  usually  neces- 
sary to  shift  the  baseband  channel  to  a more  appropriate 
location  in  the  frequency  spectrum  for  transmission.  In 
most  communications  work,  information  transmission  is  close- 
ly related  to  the  modulation  or  time  variation  of  a par- 
ticular sinusoidal  signal,  called  the  carrier,  by  a device 
called  a modulator.  In  order  to  obtain  the  information 


from  the  carrier,  the  received  signal  must  be  demodulated 
by  a similar  device  operating  as  the  inverse  of  the  modula- 
tor. Devices  used  for  modulation  and  demodulation  are 
commonly  called  MODEMS.  A typical  system  diagram  is  shown 
in  Figure  2,  where  xb(t)  is  the  baseband  signal,  x(t)  is 
the  modulated  signal,  y (t)  is  the  modulated  signal  received 
after  transmission  over  the  channel,  and  y^(t)  is  the  de- 
modulated signal  or  a replica  of  xb(t),  the  original  message 
signal.  The  modulated  signal,  y (t)  = xp(t)  + n(t)  where 
n(t)  is  noise  added  to  x (t)  by  the  channel. 


Figure  2.  Block  Diagram  of  Carrier  System 

for  Signal  Transmission  (Ref  2:119) 

Problem  and  Scope 

The  Air  Force  Communications  Service  (AFCS)  in  conjunc- 
tion with  the  Rome  Air  Development  Center  (RADC)  is  engaged 
in  research  in  the  area  of  modem  identification  for  the 
Automated  Tech  Control  (ATEC)  and  the  possible  use  of  degra- 
dation analysis  of  communications  links  based  on  the  type  of 
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modem  connected  to  the  links. 

The  goal  of  this  thesis  is  to  develop  a method  to  iden- 
tify the  type  of  modem  used  on  a communications  link  based 
on  real  time  measurement  of  y(t)  or  yfe(t)  (Figure  2).  The 
identification  consists  of  determining  the  modulation  scheme 
and  the  data  rate  at  which  the  modem  is  operating.  This 
study  is  intended  to  support  the  report  (Ref  1)  prepared  by 
AFSC/OA,  although  a theoretical  rather  than  a trial  and 
error  approach  is  taken. 

There  are  three  basic  ways  of  modulating  a sinusoidal 
carrier:  variation  of  its  amplitude,  phase,  and  frequency 

in  accordance  with  the  information  being  transmitted.  Other 
modulation  schemes  are  just  variations  and/or  combinations 
of  one  or  more  of  the  basic  modulation  schemes.  In  this 
study,  seven  different  modulation  schemes  or  modems  are 
mathematically  modeled  and  characterized  for  identification. 
They  are: 

1.  Binary  Phase  Shift  Keying  (BPSK) 

2.  Quadrature  Phase  Shift  Keying  (QPSK) 

3.  Binary  Frequency  Shift  Keying  (BFSK) 

4.  Duobinary  Frequency  Shift  Keying  (DBFSK) 

5.  Binary  Amplitude  Shift  Keying  (BASK) 

6.  Double  Sideband  Amplitude  Modulation  (DSBAM) 

7.  Double  Sideband  Amplitude  Modulation  with 

Suppressed  Carrier  (DSBAM-SC) 

The  input  signal  used  with  the  first  five  modulation 
schemes  is  a random  pulse  train  with  a basic  pulse  P(t) 
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having  magnitude  1 and  signaling  interval  T as  shown  in  Fig- 
ure 3,  while  the  input  used  with  the  last  two  modulation 


i 


schemes  is  the  special  periodic  function,  y cos  w t. 

s 


Figure  3>  Basic  Pulse  Shape 

of  Binary  Input  Signal 


The  modems  are  modeled  in  the  absence  of  noise  (i.e., 
n(t)  = 0 ) using  the  following  assumptions: 

1.  The  channel  is  ideal  and  introduces  no  anomalies 

in  the  transmitted  signal  (i.e.,  Xp  ( ) = Zp(t)  ). 

2.  The  transmitter  is  ideal  and  introduces  no  anoma- 
lies in  the  transmitted  signal. 

3.  The  input  signal  is  a random  pulse  train  repre- 

+°o 

sented  as  Xp('t)  = £ bnP(t  “ nT)  where  the  b 's 

n=o 

are  identically  independently  distributed 

Bernoulli  random  variables  each  equally  likely 

to  be  zero  or  one  and  where  the  values  of  b are 

— n 

independent  of  each  other. 

4.  The  negative  frequency  terms  do  not  overlap  the 

positive  frequency  terms  (i.e.,  X(w  - w ) does 

not  overlap  X(w  + w )). 

c 
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The  modems  are  characterized  by  the  power  spectral 
density  and  the  correlation  function  of  the  modulated  signal 


Xp(t)  shown  in  Figure  2 since  the  channel  was  assumed  ideal 
(i.e.,  ip('t)  = 2p(t)  and  yb(t)  = xb(t)  ). 


Approach  and  Outline 

First,  the  theory  and  properties  of  Fourier  transforms, 
correlation,  and  power  spectral  density  are  presented  in 
Chapter  II.  Next,  the  method  used  in  modeling  the  modems 
as  well  as  a general  operational  description  is  presented 
in  Chapter  III.  The  procedures  used  in  characterizing  the 
modems  is  presented  in  Chapter  IV.  Discussion  of  the  re- 
sults from  xhe  study  are  covered  in  Chapter  V and  the  con- 
clusions and  recommendations  are  presented  in  Chapter  VI. 


Conventions 

The  following  conventions  are  used: 

1.  The  underlining  of  a variable  indicates  that  it 
is  a random  variable  (e.g.,  x) • 

2.  The  underlining  of  an  explicit  function  of  time 
indicates  that  it  is  a random  process  (e.g., 
x(t) ) . 

3.  The  use  of  the  asterisk  (*)  as  a superscript  in- 
dicates  complex  conjugate  (i.e.,  X (w)  is  the 
complex  conjugate  of  X(w) ) . 
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II.  Fourier  Transforms . Correlation, 


and  Power  Spectral  Density- 


Introduction 

The  Fourier  transform,  correlation,  and  power  spectral 
density  find  many  applications  in  engineering  and  the  pur- 
pose of  this  section  is  to  acquaint  the  reader  with  some  of 
the  properties  of  each  of  these  used  in  this  study. 


Fourier  Transforms 

There  exists  a variety  of  notation  to  denote  the  Fourier 
transform  pair,  but  the  following  notation  and  transform  pair 
will  be  used  in  this  study.  Eq  (1)  is  the  Fourier  transform 
of  x(t),  which  in  general  is  a complex-domain  to  complex- 
domain  transformation;  whereas  Eq  (2)  is  the  inverse  Fourier 
transform  of  X(w) , which  in  general  is  a complex-domain  to 
complex-domain  transformation. 


X(w)  = F [x(t)] 

+0°  . , 

X(w)  = f x(t)e~J  Zdt  (1) 

-CO 


x(  t) 


, +°o  . 

^ f X(w)eJwtdw 


(2) 


This  transform  pair  exists  and  is  unique  provided  that  the 
function  x(t)  satisfies  the  following: 

+°o 

1.  f | x(  t ) | dt  < +oo 

— oo 

2.  Any  discontinuities  in  x(t)  are  finite  (Ref  4:9). 
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Some  Properties  of  One -Dimensional  Fourier  Transforms 

Symmetry . Magnitude  and  the  real  part  of  the  Fourier 
transform  are  even  functions,  while  phase  and  the  imaginary 
portion  of  the  Fourier  transform  are  odd  valued  functions. 
This  is  easily  seen  if  Euler's  identity  for  e~^wt  is  used, 
then  Eq  (1)  can  be  expressed  as 


X(w)  = f x(t)  (cos  wt  - j sin  wt)  dt 


where  e = cos  wt  - j sin  wt 
Now  the  real  part  of  Eq  (3)  is 


Re  [X(w)]  = f x(t)  cos  wt  dt 


while  the  imaginary  part  of  Eq  (3)  is 


Im  CX(w)]  = f x(t)  sin  wt  dt 


Because  cosine  is  an  even  function  and  sine  is  an  odd  func- 
tion, it  follows  that 


Re  [X(w)]  = Re  [X(-w)] 


Im  [X(w)]  = -Im  [X(-w)] 


Since  the  transform  is  in  general  a complex  number,  the  mag- 
nitude of  the  transform  is  even  valued  and  is  given  by 


I X(w)  | = -/{Re  [X(w)]r  + [Im  [X(w)]}' 
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The  phase  of  the  transform  is  expressed  by 

<&(w)  = tan"1  [Im  [X(w)]  / Re  [X(w)]}  (9) 

and  is,  therefore,  odd  valued  since  the  imaginary  part  of 
the  transform  is  odd. 

Linearity . The  Fourier  transform  of  a transformable 
function  is  a linear  transformation,  which  follows  directly 
from  the  additive  property  of  integration.  Let  x(t)  = 
ah1(t)  + 3h0(t)  , where  a and  3 are  constants.  Then  the 

Fourier  transform  of  x(t)  using  Eq  (1)  is 

+oO  . 

F [x(t)]  = f [och-^t)  + 0h2(t)]  e“jwtdt 

— oo 

+00  . +OQ  . 

= a f h1(t)e'Jwtdt  + 3 f h2(  t)  e"jvrcdt 

-00  -oo 

= aX1(w)  + 3X?(w)  (10) 

where  X^(w)  is  the  transform  of  h^Ct)  and  X2(^)  is  the  trans- 
form of  h2(  t)  . 

Time  Shifted.  What  happens  to  the  transform  of  a func- 
tion which  has  been  shifted  in  time?  Let  x(t)  = g(t  - T)  , 
then  the  transform  is 

f-00 

F [x(t)]  = f g(t  - T)e_i5wtdt  (11) 

-OO 

This  equation  can  be  reduced  by  using  the  change  of  variable 
u = t - T . Then  Eq  (11)  becomes 
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•f-oo 

F [x(t)]  = f g(u)e-jw(u+T)du 

— OO 

= e'*"1  J*°° g(u)e"^wudu 

— OO 

= e'^wTG(w)  (12) 

Where  G(w)  is  the  Fourier  transform  of  g(t)  by  using  Eq  (1). 
In  other  words,  the  transform  of  a function  which  has  been 
shifted  by  an  amount  T is  just  the  transform  of  the  original 
function  multiplied  by  e-JwT. 

Modulation.  What  happens  to  the  transform  of  a func- 
tion which  is  being  modulated  by  a cosinusoidal  signal? 

Let  x(t)  be  modulated  by  cos  wct,  then  the  transform  is 

+00  . . 

F [x(t)  cos  w t]  = f x(t)  cos  w te'^dt  (13) 

C -OO  c 

This  equation  can  be  reduced  by  using  Euler's  identity  for 
cos  w£t  = (e+^wc  + e-^wc)  / 2 

F [x(t)  cos  w t]  = 

V 

+ 


. Then  Eq  (13)  becomes 


+°°  jw  t _iwt 

f x(t)e  e J xdt 

-OO 

+°°  -jw  t • t 

f x(t)e  c e jwtdt 


+°°  “ j t ( w-w  ) 


1/2  f x(t)e 


dt 


+°°  -jt(w+w  ) 

1/2  J*  x(t)e  c dt  (14) 
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and  using  the  form  of  Eq  (1)  where  w = w - w and  w + w , 


then 


F [x(t)  cos  w t]  = l/2X(w  + wc)  + l/2X(w  - wq)  (15) 


Frequency  Translation.  What  happens  to  the  transform 

of  a function  which  has  been  frequency  translated?  Let 
jw  t 

z(t)  = x(t)e  c , then  the  transform  is 


F [x(t)e',Vt]  = rx(t)e'"cV>tit 


+°°  — j t ( W-W  ) 


= f x(t)e 


and  using  the  form  of  Eq  (1)  where  w = w - w , then 


jw„t 

F [x(t)e  ] = X(w  - w ) 


Multiplication . What  happens  to  the  transform  of  two 
functions  that  are  multiplied  together?  Let  x(t)  = 
x^(t)x2(t)  , then  the  transform  is 


F [xx(t)x2(t)]  = f x1(t)x2(t)e“*^Wtdt 


Using  Eq  (17)  the  result  is 


F [x1(t)x2(t)]  = f Xx(w  - wc)X2(wc)dw 


where  X-^(w  - wQ)  is  the  transform  of  x^(t)e 


and 


X2(wc)  is  the  transform  of  X2(t). 

Correlation  Function 

Definition ■ If  there  exists  a system  with  an  input  x 
and  an  output  y,  it  is  possible  to  write  a statistical  rela- 
tionship such  as  f(x,y),  the  joint  density  function  between 
x and  y.  If  there  are  no  disturbing  influences  in  the  sys- 
ten,  then  x and  y will  be  uniquely  related.  If  there  are 
disturbing  influences  in  the  system,  x and  y will  only  be 
partially  related.  The  product  moment 

+oo  +oo 

E{xy}  = f f xyf ( x, y) dxdy  (20) 

— oo  — oo 

serves  as  a measure  of  correlation  between  x and  y where 
E{xy]  is  called  the  correlation  or  expected  value  of  the 
product  of  the  random  variables  x and  y.  If  the  random 
variables  x and  y are  represented  by  random  processes 
x(t^)  and  x(t2)  respectively,  then  Eq  (20)  can  be  defined 
as 

E{x(t1)x(  t2)}  = R(t1(t2)  (21) 

where  the  function  R(t^,t2)  is  called  the  correlation  func- 
tion or  the  autocorrelation  function  of  the  random  process. 

If  a stationary  process  is  assumed,  then  instead  of 
the  ensemble  average  in  Eq  (20),  a time  average 

T 

R(x,y)  = ^ J*  x(t)y(t)dt  (22) 
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may  be  written.  Eqs  (20)  and  (22)  are  identical  if  the  pro- 
cess is  ergodic,  that  is,  one  in  which  the  ensemble  statisti- 
cal properties  do  not  vary  with  time.  All  ergodic  processes 
are  stationary,  but  not  all  stationary  processes  are  ergodic. 
Under  the  assumption  of  ergodic ity 

E{xy]  = RU.y)  (23) 

If  x(t)  and  y(t)  are  represented  by  x and  y respectively, 
then  E{x(t)y(t)}  = R[x(t)y(t)]  is  called  the  cross  correla- 
tion function  since  it  serves  as  a measure  of  the  correla- 
tion of  x(t)  and  y(t).  If  x(t)  and  y(t)  follow  each  other 
closely  and  are  of  the  same  sign  (i.e.,  they  correlate  posi- 
tively, then  E[xy}  (or  R(x,y))  will  accumulate  positively. 

If  x(t)  and  y(t)  follow  closely  and  are  of  opposite  sign 
(i.e.,  they  correlate  negatively),  then  Efxy]  (or  R(x,y)) 
will  accumulate  negatively.  On  the  other  hand,  if  x(t)  and 
y(t)  vary  independently  (i.e.,  uncorrelated)  they  tend  to 
cancel  each  other  and  the  Efxy}  (or  R(x,y))  would  approach 
zero . 

If  x(t  + t) » which  is  x(t)  delayed  by  a factor  t,  is 
substituted  for  y(t),  Eq  (22)  becomes 

R(x,x)  = i J x(t)x(t  + t) dt  (24) 

and  is  called  the  time  autocorrelation  function  and  denoted 
by  cp(r)  . 

Periodic  Functions . If  a phenomenon  which  can  be 
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designated  by  a time  history  is  repetitive  with  a constant 
duration  between  arbitrarily  selected  repetitions,  then,  it 
is  said  to  be  periodic.  The  symbolic  definition  of  a peri- 
odic function  is  expressed  as  g(t)  = g(t  + T)  . If  T is 
the  smallest  number  for  which  this  identity  holds,  it  is 
called  the  period  of  the  function.  The  relation  is  also 
assumed  to  hold  for  all  multiples  of  T.  Thus,  the  law  or 
description  of  periodic  events  is  time  invariant,  and  these 
events  can  be  classified  as  noninformative  because  they  are 
completely  predictable  without  further  observation. 

The  autocorrelation  of  a periodic  function  g(t)  is 
defined  by  the  expression 

T/2 

cp(  t)  = b f g(t)g(t  + t)  dt  (25) 

1 -T/2 


where 

T = w/2tt,  period  of  g(t) 

T = is  the  continuous  displacement  in  the  interval 
(+oo,  -oo) , independent  of  t 

A close  examination  of  the  define  •;  equation  reveals  that 
three  distinct  operations  are  involved:  displacement  by  an 
amount  r»  continuous  multiplication,  and  averaging  by  inte- 
gration over  a complete  period.  Essentially,  the  displace- 
ment and  averaging  provide  a means  of  transforming  the  func- 
tion to  a new  time  domain.  This  process  also  eliminates 
only  phase  information,  and  only  the  magnitude  information 
can  then  be  transformed  to  the  frequency  domain. 
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(26) 


The  Fourier  expansion  of  g(t)  is 

+oo  -jw  t 

g(t)  = S G(n)e 

n=-<» 


where 


G(n)  = ^ - jbn 


T/2 

s 

-T/2 


-jnwt 

g(t)e  n dt 


(27) 


The  Fourier  expansion  of  g( t + t)  Is 


+°°  -jnw  (t+j) 

g(t  + t)  = 2 G(n)e 

n= 

Substituting  Eq  (28)  into  Eq  (25) 

T/2  -i  +°°  jnw  (t+r) 

cp( t)  = b f g(t)  b 2 G(n)e  n dt 
1 -T/2  n=-°° 


(28) 


(29) 


or  rearranging  Eq  (29) 


, +oo  jnw  r i T/2  t 

cp(T)  2 G(n)e  £ S g(t)e 

n=-»  -T/2 


dt  (30) 


The  integral  in  Eq  (30)  is  the  complex  conjugate  of  Eq  (2?). 
Substituting  G*(n)  into  Eq  (3) 


,_X  1 r CG(n)G*(n)]  eJnWnT 

<P(T)  T T 


~ ' 0(n)  I 2 BJnV 


— £ 

T n--»  T 


(3D 
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where  | G(n) | = G(n)G  (n)  . The  autocorrelation  function 

with  an  argument  of  zero  (r  = 0)  is 

T/2 

Cp(  0)  = i f g2(t)dt  (32) 

1 -T/2 

which  is  the  mean  square  value  of  g(t).  If  it  is  assumed 
that  there  is  a load  resistor  of  one  ohm  and  that  g(t)  rep- 
resents either  current  or  voltage,  then  the  mean  square 
power  consumed  is  the  sum  of  the  power  contributed  by  the 
individual  harmonics  into  which  g(t)  has  been  resolved. 

This  is  known  as  Parseval's  Theorem  for  periodic  functions. 

If  |G(n)|  / T is  defined  as  the  power  spectrum  $(n)  , 

then  from  Eq  (31) 

+oo  jnw  t 

cp(  T ) = if  Z $(n)e  n (33) 

n=-°° 


and  inversely 


$(n) 


T/2 

S, 

-T/2 


cp( 


-jnw  T 

t)  e 


= | G(n) | 2 / T 


(34) 


Eqs  (33)  and  (34)  form  the  autocorrelation  theorem  for  a 
periodic  function.  That  is,  the  autocorrelation  function 
and  the  power  spectrum  are  Fourier  transforms  of  each  other. 

In  summary,  some  important  properties  covered  in  this 
section  on  periodic  functions  (Ref  6:1166)  are: 

1.  A periodic  function  can  be  represented  as  the 
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summation  of  discrete,  ordered  harmonics  having 
the  properties  of  amplitude  and  phase  relations 
in  the  time  domain. 

2.  These  properties  can  be  transformed  to  the  fre- 
quency domain. 

3.  Autocorrelation  transforms  only  the  magnitude 
properties  of  the  unordered  harmonics  to  a new 
time  domain. 

4.  These  magnitude  properties  can  be  transformed  to 
the  frequency  domain. 

Random  Functions . If  a phenomenon  which  can  be  desig- 
nated by  a time  history  is  repetitive  with  varying  duration 
between  arbitrarily  selected  repetitions,  then  it  is  said  to 
be  random.  Thus,  the  law  or  description  of  random  events  is 
time  variant  or  unordered  in  time  . These  events  are  classi- 
fied as  informative  because  they  are  predictable  only  in  a 
statistical  sense.  Unlike  the  periodic  phenomenon  which  is 
completely  predictable  from  the  assumed  exact  function  g(t), 
the  random  phenomenon  can  only  be  predicted  from  an  assumed 
exact  probability  density  function.  The  exactness  assump- 
tion for  the  periodic  function  is  usually  obvious,  although 
in  the  end  it  is  dependent  on  induction.  The  exactness 
assumption  for  probability  densities  is  not  so  obvious,  and 
it  is  dependent  on  such  concepts  as  stationarity  and  ergod- 
icity  as  well  as  induction.  However,  these  difficulties  are 
easily  removed  by  defining  the  autocorrelation  function  for 
the  random  function  as  a limiting  process  for  an  infinite 
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observation.  Thus,  the  observation  must  be  long  enough  to 
establish  the  exactness  of  the  probability  densities  in- 
volved in  correlation,  or  for  lesser  lengths  it  can  be  as- 
sumed that  an  approximation  is  made  of  the  exact  autocorrela- 
tion function. 

With  the  above  concepts,  the  random  autocorrelation 
function  can  be  derived  directly  from  the  periodic  autocor- 
relation function  by  allowing  the  period  to  go  to  infinity. 
From  Eq  ( 25) 


T/2 

cp(r)  = f g(t)g(t  + T)dt  (35) 

1 -T/2 


As  T-»+°°  and  cp— rp-^ , the  random  autocorrelation  is  given  as 


T 

CPj/t)  = 2T  £ Si(t)g;L(t  + r)dt  (36) 


where  g^(t)  is  considered  a stationary  random  function. 

Since  the  random  function  g^(t)  cannot  have  a periodic  com- 
ponent by  definition,  it  can  be  shown  that  the  random  auto- 
correlation function  must  decay  to  zero  and  is,  therefore, 

+oo 

aperiodic  and  finite.  Now  the  / |cp-.(T)|dT  is  also  fi- 

-oo 

nite  and  this  satisfies  the  bounding  assumption  for  a Fourier 
analysis . 

The  representation  of  a random  autocorrelation  function 
in  the  frequency  domain  can  also  be  derived  from  the  peri- 
odic approach.  From  Eqs  (33)  and  (3*0 
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cp(r)  = E e 


+oo  jnw  r , T/2  -jnTW 

v-  ~ XI  X « _ _ / _ \ _ I 


f cp(r)e 
1 -T/2 


As  T-+«f  l/T-dw/2-rT , nwn-w  and  cp(t)-'P1(t)  then 


CD^r)  = 2n  / eJWTdw  f cp1(T)e" JWTdr 


As  before,  the  now  continuous  spectrum  of  cp-j/t)  is  defined 


$l(w)  = f cp1(T)e-JWTdT 


and  it  forms  a Fourier  transform  pair,  in  integral  form, 


with 


cp-]_(t)  = 2n  f $1(w)e^WTdw  (40) 

-oo 

The  physical  significance  of  ^(w)  can  be  shown  by  having 
T = 0 in  Eq  (40) . Then 

1 +°° 

Cp^O)  = 2tt  J*  ^i(w)dw  (41) 

-OO 

From  the  defining  equation  for  the  autocorrelation  function, 
Eq  ( 36) , when  r = 0, 


<Pl(°)  = IT  £ CSq(  t)]^dt 
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which  is  the  mean  square  value  of  the  function  g^(t) . 
Therefore,  from  Eqs  (41)  and  (42) 


T o ■,  +oo 

It  f Cg!(t)]2dt  = J*  $1(w)dw 
-T  — oo 


(43) 


From  Eq  (43),  $-^(w)  is  referred  to  as  the  power  density 
spectrum  or  power  spectral  density  of  the  random  function 
g-j_(t).  If  gj_(t)  represents  a voltage  across,  or  a current 
through,  a 1 ohm  resistor,  then  the  mean  value  of  g-^(t)  is 
the  mean  power  taken  by  the  load.  This  mean  power  consumed 
is  equal  to  the  integral  of  <J>^(w)  with  respect  to  the  angu- 
lar frequency  w over  the  entire  range  of  frequencies.  Con- 
sequently, $-^(w)  represents  the  power  spectral  density  of 
g-^(t)  (Ref  7:58).  Since  (^(w)  is  defined  as  the  Fourier 
transform  of  cp-^(r)  , it  is  a continuous  spectrum  of  the  auto- 
correlation function. 

Some  important  properties  of  random  functions  are: 

1.  A random  function  cannot  be  represented  as  the 
summation  of  discrete  harmonics,  but  it  can  be 
statistically  represented  as  the  integrated  effect 
of  a continuous  band  of  unordered  harmonics. 

2.  Autocorrelation  transforms  the  amplitude  proper- 
ties of  these  harmonic  bands  to  a new  finite 
time  domain. 

3.  These  amplitude  properties  can  be  transformed  to 
the  frequency  domain. 
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Power  Spectral  Density 

Periodic  Functions . The  nomenclature,  power  spectral 
density,  does  not  actually  exist  for  periodic  functions. 

The  correct  nomenclature  for  periodic  functions  is  power 
spectrum  and  is  represented  hy  Eq  (34).  In  this  equation, 

p 

| G-j^C n)  | is  the  square  of  the  magnitude  spectrum  of  g(t), 
the  periodic  function. 

Random  Functions . The  power  spectral  density  for  a 
random  function  is  given  by  Eq  (39).  where  rp^(r)  is  the 
autocorrelation  function  of  a random  function.  The  auto- 
correlation function  is  represented  by  Eq  (40)  and  with  Eq 

(39)  forms  a Fourier  transform  pair.  However,  Eqs  (39)  and 

(40)  hold  only  when  the  random  function  is  stationary  in  a 
wide  sense  (i.e.,  the  random  function's  statistical  proper- 
ties are  invariant  to  time  translation) . A more  general 
form  for  the  power  spectral  density  is  given  as 

$x(w)  = E{|G(w)| 2} 

= E{|G(w)G*(w)| } (44) 


where 
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Modem  Model  Description  and 
Theory  of  Operation 


Before  an  information-bearing  signal  is  transmitted 
through  a communications  channel,  some  type  of  modulation 
process  is  typically  utilized  to  produce  a signal  which  can 
easily  be  accommodated  by  the  channel. 

The  logical  choice  of  modulation  techniques  is  influ- 
enced by  the  characteristics  of  the  message  signal,  the  char- 
acteristics of  the  channel,  the  performance  desired  from  the 
overall  communication  system,  the  use  to  be  made  of  the 
transmitted  data,  and  the  economic  factors  which  are  always 
important  in  practical  applications.  For  this  study,  the 
modulation  techniques  chosen  were  influenced  by  their  ease 
of  modeling  and  use  in  present  day  communications  systems. 

There  are  two  basic  types  of  modulation:  binary  and 
analog.  In  the  binary  case  there  are  three  basic  ways  of 
modulating  a sinusoidal  carrier:  variation  of  its  amplitude, 
phase,  and  frequency  in  accordance  with  the  information 
being  transmitted.  Basically,  this  corresponds  to  switch- 
ing these  three  parameters  between  either  of  two  possible 
values.  The  three  basic  types  of  binary  modulation  are: 
amplitude  shift  keying  (ASK),  phase  shift  keying  (PSK),  and 
frequency  shift  keying  (FSK) . 

In  the  analog  case  there  are  two  types  of  modulation: 
continuous-wave  and  pulse.  In  continuous-wave  modulation, 
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one  of  the  three  parameters  is  again  varied  proportionally 
to  the  message  signal  such  that  a one-to-one  correspondence 
exists  between  the  parameter  and  the  message  signal.  The 
carrier  is  usually  assumed  sinusoidal,  though  not  necessar- 
ily so.  The  three  basic  types  of  continuous  wave  modulation 
are:  amplitude  modulation  (AM),  phase  modulation  (PM),  and 

frequency  modulation  (FM) . 

In  analog  pulse  modulation,  some  characteristic  of  a 
pulse  is  made  to  vary  in  one-to-one  correspondence  with  the 
message  signal.  Since  the  pulse  is  characterized  by  three 
parameters  - amplitude,  width,  and  position  - there  are 
three  types  of  analog  pulse  modulation.  They  are:  pulse- 
amplitude  modulation  (PAM),  pulse-width  modulation  (PWM), 
and  pulse-position  modulation  (PPM) . 

Other  modulation  techniques  exist  and  are  variations 
and/or  combinations  of  the  modulation  techniques  previously 
listed.  Different  modulation  techniques  can  also  be  real- 
ized by  changing  the  carrier  signal  from  sinusoidal  to  one 
that  is  pulsed. 

In  this  study  the  carrier  is  assumed  sinusoidal  and 
five  cases  of  binary  modulation  and  two  cases  of  analog 
modulation  are  considered.  The  different  modulation  tech- 
niques or  modems  modeled  are: 

1 . BPSX 

2 . QPSK 

3.  BFSK 

4 . DBFSK 
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5 . BASK 

6 . DSBAM 

7.  DSBAM -SC 


BPSK  Modem 

In  BPSK  the  phase  of  the  carrier  is  switched  between 
two  values  according  to  the  two  levels  of  the  digital  input 
signal  to  the  modem.  The  two  values  are  usually  separated 
by  tt  radians  corresponding  to  the  two  signal  levels.  (This 
is  the  only  choice  considered  in  this  report.)  This  tech- 
nique is  sometimes  called  phase  reversal  keying  (PRK) . 

The  effect  of  phase  reversal  modulation  is  to  produce 
a double  sideband  suppressed-carrier  amplitude -modulated 
waveform  where  the  information  signal  is  digital.  Thus,  the 
PSK  signal  can  be  written  as 

xi(t)  = Acxd(t)  cos  wct  i = 0,1  (45) 

bn 

where  xd(t)  = (-1)  and  has  value  +1  or  -1  according  to 
the  input  digital  code  (b  = 0 or  1)  during  any  bit  inter- 
val having  a duration  T,  Ac  is  the  peak  carrier  amplitude, 
and  wc  is  the  carrier  frequency.  If  the  carrier  frequency, 

w , is  an  exact  multiple  of  the  bit  rate,  then  w = nBu/T  , 
c c 

where  the  bit  rate  is  T-1.  The  waveforms  of  the  modulated 
signal,  xQ(t)  and  x-^(t),  can  be  written  as 


X0(t)  = Accos  wct 

(46) 

xx(t)  = -xQ(t) 

(47) 
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A simple  model  of  a BPSK  modem  is  shown  in  Figure  4. 
The  switch  in  Figure  4 operates  at  a rate  determined  by  the 
signaling  interval,  T,  of  the  binary  input  signal. 
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Figure  4.  Model  for  a BPSK  Modem 


QPSK  Modem 

With  binary  digital  communications  systems,  one  of  only 
two  possible  signals  can  be  transmitted  during  each  T-second 
signaling  interval.  In  an  M-ary  system,  one  of  M possible 
signals  may  be  transmitted  during  each  T-second  signaling 
interval,  where  M > 2.  Thus,  binary  data  transmission  is 
a special  case  of  M-ary  data  transmission.  Each  possible 
transmitted  signal  of  an  M-ary  message  sequence  is  referred 
to  as  a character  or  symbol.  The  rate  at  which  M-ary  sym- 
bols are  transmitted  through  the  channel  is  called  baud 
rate  in  bauds. 

One  of  the  most  common  M-ary  systems  is  four-phase,  or 
quadriphase  shift  keying  (QPSK). 

In  general  a QPSK  modem  can  be  thought  of  as  two  BPSK 
modems  in  parallel  in  which  the  separate  carriers  are  in 
phase  quadrature.  Thus,  for  a given  baud  rate,  T-1,  in 
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bauds  the  rate  at  which  binary  digits  are  sent  through  the 
channel  for  QPSK  is  double  the  rate  that  can  be  achieved 
with  a binary  modem  operating  at  the  same  signaling  rate. 
Thus,  the  signal  from  the  QPSK  modem  has  a signaling  inter- 
val twice  that  of  the  binary  input  signal.  The  QPSK  signal 
can  be  written  as 

xi(t)  = xq(t)  + Xp(t)  i = 0,1, 2,3  (48) 

where 

Xq(t)  = Vdq(t)  Sin  Wct  (49) 

and 

x.p(t)  = Acxdp(t)  cos  wct  (50) 

Eq  (^9)  represents  the  quadrature  phase  modulated  signal  and 
the  in-phase  modulated  signal  is  represented  by  Eq  (50). 

The  carriers,  cos  w t and  sin  wct  are  modulated  by  the  bit 
streams  xdq(t)  and  xdp(t)  which  are  obtained  by  grouping  the 

bits  of  the  binary  signal  xq(t)  with  half  the  bit  period  of 

xdq  and  xdp,  two  bits  at  a time.  The  values  for  xd(t)  are 
-1  or  +1  and  are  determined  by  the  input  digital  code,  b . 
Therefore,  from  the  above  discussion,  xdq(t)  and  xdp(t)  have 
values  +1  or  -1  depending  on  the  digital  codes  bn  and  bm,  re- 
spectively. The  bn  and  bm  digital  codes  result  from  the  cod- 
ing of  bQ.  Thus,  xdq(t)  = (-1)  n and  xdp(t)  = (-1)  m 

with  bn  and  bm  each  having  value  of  0 and  1.  Since  the 
phase  of  the  transmitted  signal  is  9^  = tan-1  ( xdq  / xdp) 
it  can  be  seen  that  9^  can  take  on  four  possible  values, 

±45°  and  ±135° • Thus,  four  signals  can  be  obtained  as  shown 
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by  the  substitution  of  Eqs  (49)  and  (50)  into  Eq  (48)  when 
xd^(t)  and  xdp(t)  each  have  values  of  +1  or  -1. 

A simple  model  of  a QPSK  modem  is  shown  in  Figure  5- 


Figure  5-  Model  for  a QPSK  Modem 

BFSK  Modem 

Generating  a carrier  waveform  by  BFSK  consists  of 
switching  between  two  constant  frequency  sources  with  the 
source  selected  being  dependent  on  the  binary  input  signal 
which  is  coded  with  zeros  and  ones.  Thus,  the  transmission 
of  a code  zero  may  correspond  to  a pulse  of  frequency 
w + aw  , while  a code  one  selects  a frequency  w„  - aw  , 
where  w is  the  carrier  frequency  and  aw  is  the  peak  fre- 
quency  deviation.  Thus,  the  FSK  signal  can  be  written  as 

xi(t)  = Ac  cos  (wc  + xd(t)AW)t  1 = 0.1  (5D 
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where  xd(t)  = (-1)  and  has  value  of  +1  or  -1  according 
to  the  digital  code  (b  = 0 or  1)  during  any  bit  interval 
having  a duration  T,  and  AQ  is  the  peak  carrier  amplitude. 
The  waveforms  xQ(t)  and  x1(t)  can  be  written  as 

x0(t)  = Ac  cos  (wc  + aw) t 

= Ac  cos  wxt  (52) 


x-^t)  = Ac  cos  (wc  - aw)  t 

= Ac  ccs  w2t  (53) 

where 

w,  = w + aw 
1 c 

W2  = Wc  - AW 

A simple  model  of  a BFSK  modem  is  shown  in  Figure  6. 

The  switch  operates  at  a rate  determined  by  the  signaling 
interval,  T,  of  the  binary  input  signal. 


Figure  6.  Model  for  a BFSK  Modem 
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DBFSK  Modem 


i 


In  generating  a DBFSK  waveform,  the  binary  input  data 
is  encoded  (duobinary  encoding)  before  being  modulated. 
Figure  7 shows  one  of  the  ways  in  which  the  binary  input 
signal  is  encoded. 


Figure  7.  Duobinary  Encoder  (Ref  5:400) 


Binary  data  pulses  are  applied  to  one  input  of  a sim- 
ple AND  gate,  and  a stream  of  binary  "clock"  pulses  of  the 
same  repetition  rate  as  the  data  pulses  is  applied  to  the 
other  input.  Whenever  a data  space  (zero  bit)  occurs,  the 
AND  gate  yields  an  output  pulse  which  is  applied  to  a "flip- 
flop"  circuit  or  bi-stable  multivibrator.  The  multivibrator 
has  the  characteristic  of  changing  its  output  voltage  from 
zero  to  plus  or  plus  to  zero,  whenever  it  is  triggered  by 
the  AND  gate.  Thus,  for  every  data  space,  the  output  of 
the  flip-flop  changes  state.  Whenever  a data  mark  is  re- 
ceived, no  change  occurs.  The  binary  waveform  from  the 
flip-flop  is  then  passed  through  a low-pass  shaping  filter. 

J 
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Because  of  the  encoding  process,  the  signal  now  has  one  very 
unique  property.  As  it  passes  through  the  shaping  filter, 
the  waveform  is  altered  so  that  it  tends  to  occupy  three 
amplitude  levels.  Superficially,  this  three-level  signal 
resembles  a bipolar  signal  as  shown  in  Figure  8.  This  duo- 
binary signal  has  one  unique  feature  in  that  the  signal  can- 
not transition  from  one  extreme  level  to  the  opposite  ex- 
treme within  one  bit  interval. 


Figure  8.  Resulting  Waveforms  of  a Binary 
Signal  Inserted  into  the  Diagram 
Shown  in  Figure  7 (Ref  5:^00) 


A continuous  series  of  input  data  spaces  causes  the 
flip-flop  to  change  state  repeatedly,  thus  creating  a series 
of  binary  pulses  at  the  data  rate.  Since  this  is  too  fast 
to  pass  through  the  shaping  filter,  a zero  level  results. 

A continuous  series  of  marks  results  in  the  flip-flop  re- 
maining in  whichever  state  it  was  in  at  the  beginning  of 
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the  sequence.  The  single  voltage,  tone,  or  phase  by  which 
this  marking  condition  is  transmitted  will  pass  through  the 
filter.  Since  either  plus  or  minus  represents  a mark,  the 
signal  will  not  change  for  consecutive  marks,  but  only  when 
a space  follows  a mark.  A second  space  holds  the  output  in 
the  zero  condition,  but  a mark  following  the  space  returns 
the  signal  to  a mark  condition  as  shown  in  Figure  8. 

The  duobinary  transformation  (Ref  11)  used  in  this 
study  is 

H(w)  = Jl/2  cos  (wT/2)  (54) 

where  T-"'"  is  the  data  rate  and  H(w)  is  a transformation  in 
the  frequency-domain.  Therefore,  the  duobinary  signal  of 
Figure  7 can  be  represented  as 

z(t)  = h(t)xa(t)  (55) 

where 

z(t)  = is  the  duobinary  signal 

h(t)  = is  the  inverse  Fourier  transform  of  H(w) 

x_(t)  = is  the  binary  input  signal 

a 

After  the  input  waveform  has  been  encoded,  the  result- 
ing duobinary  signal  is  modulated  using  the  FSK  modulation 
scheme.  However,  instead  of  switching  between  two  constant 
frequency  sources,  the  switching  occurs  between  three  con- 
stant frequency  sources  with  the  source  selected  being  de- 
pendent on  the  duobinary  signal.  Since  the  duobinary  sig- 
nal has  three  levels  (in  this  case  represented  by  +1,  0, 
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and  -1) , the  transmission  of  a code  +1  corresponds  to  a 

pulse  of  frequency  of  wc  + aw.  For  a code  0 the  frequency 

is  w and  for  a code  -1  the  frequency  is  w - aw.  Thus,  the 
c c 

DBFSK  signal  can  be  written  as 

xi(t)  = Ac  cos  (wc  + xdb(t)Aw)t  i = 0,1,2  (56) 

where  xdb(t)  has  values  +1,  0,  and  -1,  depending  on  the  duo- 
binary coding  of  the  binary  input  signal.  The  waveforms 
then  are 


xQ(t)  = Ac  cos  (wc  + Aw)t 

= Ac  cos  w1t  (57) 

Xf(t)  = Ac  cos  wct 

= Aq  cos  w2t  (58) 

x2(t)  = Ac  cos  (wc  - Aw ) t 

= Ac  cos  w3t  (59) 


where 


W1  = 

wc  + 

AW 

when 

xdb(t) 

is 

+ 1 

II 

CM 

s 

wc 

when 

xdb(t) 

is 

0 

w3  “ 

w + 
c 

AW 

when 

xdb(t) 

is 

-1 

A simple  model  of  a DPFSK  modem  is  shown  in  Figure  9. 
The  switch  operates  at  a rate  determined  by  the  signaling 
interval,  T,  of  the  duobinary  input  signal. 
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Figure  9.  Model  for  a DBFSK  Modem 


BASK  Modem 

In  BASK  the  amplitude  of  a carrier  is  switched  between 
two  levels,  usually  the  extremes  of  full  on  and  full  off. 

As  a result,  ASK  is  sometimes  referred  to  as  on-off  keying 
(00K).  The  on  condition  might  typically  correspond  to  a 
code  one,  while  off  corresponds  to  a code  zero.  Thus,  the 
ASK  signal  is  written  as 

A 

xi(t)  = [1  + xd(t)]  cos  wct  i = 0,1  (60) 

bn 

where  xd(t)  = ( -1;  and  has  value  +1  or  -1  according  to 
the  digital  code  (b  = 0 or  1)  during  any  bit  interval 
having  a duration,  T,  and  A is  the  peak  carrier  amplitude. 
The  "1"  in  Eq  ( 60)  is  the  dc  level  added  to  the  digital  sig- 
nal xd(t).  The  waveforms  xQ(t)  and  x-^t)  can  be  written  as 

xQ(t)  = Ac  cos  wct  (61) 
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x1(t)  = 0 (62) 

A simple  model  of  a BASK  modem  is  shown  in  Figure  10. 
The  switch  operates  at  a rate  determined  by  the  signaling 
interval,  T,  of  the  binary  input  signal. 


k 

1 

2 

Cos 

^Li) 

« 

r^j 

o 

Figure  10.  Model  for  a BASK  Modem 


DSBAM  Modem 

In  DSBAM  the  amplitude  of  the  carrier  is  varied  propor- 
tionally to  the  message  signal  such  that  a one-to-one  cor- 
respondence exists  between  the  carrier  amplitude  and  the 
message  signal.  Since  DSBAM  results  when  a dc  bias  is 
added  to  the  message  signal  prior  to  the  modulation  process, 
the  resulting  waveform  is 

x(t)  = Ac  [A  + xm(t)]  cos  wct  (63) 

where 

A = maximum  amplitude  of  the  carrier 
wc  = carrier  frequency 
xJn(t)  = analog  message  signal 
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A = dc  bias  added  to  message  signal 
A simple  model  of  a DSBAM  modem  is  shown  in  Figure  11. 


Figure  11.  Model  for  a DSBAM  Modem 


DSBAM-SC  Modem 

The  only  difference  between  DSBAM  and  DSBAM-SC  is  that 
in  Eq  (63),  the  dc  bias,  A,  is  zero.  Thus,  the  resulting 
DSBAM-SC  waveform  is  represented  as 

x(t)  = Acxm(t)  cos  wct  (64 

A simple  model  of  a DSBAM-SC  modem  is  shown  in  Figure 


IV . Analytical  Procedures  for 
Modem  Characterization 

Introduction 

In  this  chapter  the  analytical  procedures  for  character- 
izing the  modems  modeled  in  Chapter  III  are  formulated.  Each 
modem  is  characterized  primarily  by  determining  the  power 
spectral  density  and  secondarily  by  determining  the  autocor- 
relation function  (except  for  the  DBFSK,  DSBAM,  and  DSBAM-SC 
modems)  of  the  output  signal,  xp(t),  shown  in  Figure  2. 

The  input  signal  to  each  of  the  binary  modems  is  a ran- 
dom process  and  for  the  binary  modems  this  input  has  only 
two  possible  values,  namely,  zero  and  one.  The  binary  sig- 
nal is  a random  pulse  train  with  each  pulse  having  period  T . 

A sample  pulse  is  shown  in  Figure  3-  The  output  signal  from 
each  of  the  binary  modems  modeled  is  also  a random  process. 

The  input  signal  to  the  analog  modems  is  a special  per- 
iodic function  given  by  (Ref  12:428) 

a(t)  = v cos  wgt  (65) 

where 

w_  = frequency  of  the  special  function  and  w < w„ 

Y = amplitude  factor  0 < y < 1 

The  expected  or  mean  value  of  the  output  signal  will 
also  be  determined  (except  for  DBFSK,  DSBAM  and  DS3AM-SC 
modems)  as  it  will  be  needed  in  determining  the  autocorrela- 
tion function.  Once  the  autocorrelation  function  is 
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calculated  then  it  can  be  determined  whether  the  output  sig- 
nal of  the  modem  being  analyzed  is  a nonstationary  or  sta- 
tionary process  in  the  wide  sense.  A wide  sense  stationary 
process  is  one  that  has  a mean  that  is  constant  or  zero,  and 
whose  second  moments  are  invariant  to  time  translation. 

This  invariance  implies  that  the  underlying  physical  mecha- 
nism producing  the  process  in  not  changing  with  time.  The 
determination  of  whether  the  output  signal  is  a stationary 
process  in  the  wide  sense  or  a nonstationary  process  has  a 
bearing  on  how  the  autocorrelation  is  related  to  the  power 
spectral  density.  If  the  process  is  not  stationary  then 
Eq  (44)  is  used  instead  of  Eq  (39)  to  find  the  power  spec- 
tral density. 

BPSK  Modem  Characterization 

Using  the  modem  model  shown  in  Figure  4 of  Chapter  III 
with  an  input  a random  pulse  train,  the  output  signal  of  a 
BPSK  modem  is  represented  as 


N b 

x (t)  = £ (-1)  cos  w tP(t  - nT)  (66) 

p n=  o 


where 

random  variable  equally  likely  to  be  one 
or  zero 

the  n signaling  interval  of  the  input 
signal  summed  over  n intervals  with  each 
pulse  having  period  T 

wc  = carrier  frequency 


N 

£ P(t  - nT)  = 
n=o 
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The  mean  or  expected  value  of  Xp(t)  is 

, N . 

E{ x (t)}  = E £ (-1)  n cos  w tP(t  - nT)  (67) 

p ' n=o  c 1 

Now  using  the  definition  of  expectation  (Ref  10:138-140),  Eq 
(67)  can  be  expressed  as 

N b 

E(Xp(t)}  = £ cos  wctP(t  - nTjEfC-l)-11}  (68) 

since  the  expectation  of  a constant  times  a random  variable 
is  the  constant  times  the  expectation  of  the  random  variable. 
Therefore,  using  the  given  assumption  about  the  input  signal 
having  zeros  and  ones  equally  likely,  the  probability  P of  a 
one  or  zero  occurring  is  P(bn  = 1)  = P(bn  = 0)  = 1/2 
Then,  substituting  in  for  b^,  the  values  zero  and  one,  and 
using  the  definition  of  expectation  of  a sum  of  discrete 
random  variables  (Ref  10:138) 

E((-irn]  = P(bn  = 0)(-l)°  + P(bn  = 1)  ( -l)1 
= ( l/2 ) ( 1 ) + ( 1/2) ( -1) 

= 0 (69) 

Therefore,  substituting  the  result  of  Eq  (69)  into  Eq  (68) 

E[xp(t)}  =0  (70) 

results  in  the  expected  value  or  mean  of  the  output  signal 
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to  be  zero 

The  autocorrelation  function  of  Xp(t)  using  Eq  (21)  is 
Rttp.tg)  = E{xp( t1)xp(t2)} 

N N 

= E E E cos  w t,  cos  w t pP ( t-i  - nT) 
n=o  m=o  c x c ^ 

b b 

P(t2  - mT)(-l)  n(-l)  m (71) 

and  using  the  definition  of  expectation,  Eq  (70)  is  ex- 
pressed as 

N N 

R(tp,t0)  = E E cos  wct^  cos  wct2P( t1  - nT) 
n=o  m=o 

b b 

P(t2  - mT)E{ ( -1)  n(-l)  m}  (72) 

In  Eq  (72),  the  output  signal,  Xp(t),  is  sampled  at  two  dif- 
ferent times,  t1  and  t2»  which  may  or  may  not  lie  within  the 

same  signaling  interval,  T.  The  signaling  intervals  are  de- 

b^  bm 

termined  by  the  values  of  m and  n.  Now  the  E{ ( -l)-n( -1)—  } 
in  Eq  (72)  is  found  by  referring  to  Tables  I and  II. 


TABLE  I 

Comparison  of  Sample  Pulse  Values  and  the 
Probability  of  the  Pulses  Occurring  (n/m) 
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TABLE  I 


Comparison  of  Sample  Pulse  Values  and  the 
Probability  of  the  Pulses  Occurring  (n=m) 


Sample  Pulse  Value 

Probability 

of  Occurrence 

— n 

— m 

P<V 

P(bm) 

' — m7 

0 

0 

1/2 

1/2 

1 

1 

1/2 

1/2 

These  two  tables  represent  the  two  cases  for  determining 
b b 

the  E{(-1)  (-1)  },  namely,  when  the  sample  times,  t-^  and 

t2,  occur  in  different  signaling  intervals,  n / m,  and  when 
the  sample  times,  t^  and  t2,  occur  in  the  same  signaling 
interval,  n = m.  Using  the  definition  for  expectation  of 
a sum  of  discrete  random  variables  and  remembering  that  the 
pulses  are  independent,  then  for  the  case  when  n / m 

E[(-l)"n(-l)“m}  = E{(-irn}E{(-irm} 

= P^  = 0)P(bm  = ))(-l)°(-l)° 

+ P(bn  = 0)P(bm  = 1)(-1)°(-1)1 

+ P(bn  = l)P(bm  = 0)(-l)1(-l)° 

+ P(bn  = l)P(bm  = 1)(-1)1(-1)1 

= 1/4  - 1/4  - 1/4  + 1/4 

= 0 (73) 

In  the  case  where  n = m,  only  two  conditions  occur,  namely 
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when  n = id  = 0 and  n = m = 1.  As  a result 

E{(-l)"n(-l)“m}  = E{(-l)"n}E{(-l)"m] 

= PC^  = 0)(-l)°(-l)° 

+ P(bn  = 1)(-1)1(-1)1 

= 1/2  + 1/2 

=1  (74) 

Therefore,  substituting  the  result  of  Eq  (73)  into  Eq  (72), 
the  autocorrelation  function 

H(t1,t2)  = o (75) 

when  n / m or  when  t1  and  t2  are  not  in  the  same  signaling 
interval.  Substituting  the  result  of  Eq  (74)  into  Eq  (72) 

N N 

R(tlft2)  = S E cos  wQt1  cos  wct2P(t1  - nT)P(t2  - mT)  (76) 
n=o  m=o 

when  n = m or  when  t^  and  t2  are  in  the  same  signaling  in- 
terval. Since  n = m in  Eq  (76),  it  can  be  reduced  to 

R(t1,t2)  = cos  wctl  cos  wct2  (77) 

N N 

as  the  term,  E E P(t,  - nT)P(t2  - nT)  = 1 
n=o  m=o 

Combining  Eqs  (75)  and  (77) 
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R( tx,t2) 


0 when  tp  and  t2  are  not  in 

the  same  signaling  interval 

cos  v^t^  cos  wct2  when  t^  and  t2  are  in  the 

same  signaling  interval 


(78) 


From  the  definition  of  a wide-sense  stationary  process,  it 
can  be  seen  that  x.(t)  is  nonstationary  as  it  is  time  de- 

v 

pendent.  Therefore,  Eq  (44)  is  used  to  find  the  power  spec- 
tral density  of  the  output  signal  from  the  modem. 

In  order  to  find  the  power  spectral  density,  the  output 
signal  Xp(i;)  must  be  transformed  to  the  frequency-domain. 
Using  the  properties  of  the  Fourier  transform  outlined  in 
Chapter  II,  the  frequency-domain  transform  of  Xp(t)  is  rep- 
resented as 


Xp(w)  = F [Xp(t)] 


= F 


N b 

E (-1)  cos  w tP(t  - nT) 
n=o 


(79) 


where 


F [cos  wQt]  = 1/26 (w  + wc)  + 1/26 (w  - wc) 


and 


F [P(t  - nT)] 


T sin  (wT/2 )e-3»T(n+l/2) 
(WT/2) 


(80) 


(81) 


Therefore, 


41 


Vw) 


N b 

= E (-1)  n(T/2) 


sin  [(w  + w )T/2]  e 
(w  + wc)T/2 


■ jT(w+w  ) (n+l/2) 


- jT(w-w  ) (n+l/2) 

sin  [(w  - wc)T/2]  e 
(w  - w )t/2~ 


Now  from  Eq  (44) 


|Xp(w)l 2 = Xp(w)X*(w) 


E[|Xp(w)|2}  = E{Xp(w)Xp(w)} 


where  X*(w)  is  the  complex  conjugate  of  Xp(w)  . Since  xp(t) 
was  assumed  real,  then  Xp(w)  i-s  just  the  conjugate  of  the 
real  process  X (w) • The  power  spectral  density  is  then 


o N N b b 9 

EClX^Cw)!2}  = E E E (-1)  n(-D  m(T2/4) 
p n=  o m=  o 


sin  [(w  + w )T/2]  e 

~ -w  

(w  + wc)t/2 


-jT(w+w  ) (n+1/2) 


sin  [(w  + w )t/2]  e 
(w  + wc)T/2 


sin  [(w  - wc)T/2]  e 
(w  - wc)t/2 


- jT(w+w  ) (m+l/2) 

sin  [(w  - wc)t/2] 
+ (w  - wc)T/2 


-jT(w-w  ) (n+1/2)  -jT(w-w  ) (m+ 1/2)1 
c e c 
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or  combining  terms 


~ NN, 
E{|X  (w)H  = E E T /4 
p n=o  m=o 


o 

sin  [(w  + w )t/2]  e 
C(w  + wjT/2]2 


-jT(w+w  ) (n+l/2)  -jT(w+w  ) (m+l/2) 


sin^  [(w  - wc)T/2]  e 
[(w  - w )T/2] 2 


-jT(w-w) (n+l/2)  -jT(w-w) (m+l/2) 


E{(-l)bn(-l)bnl}  (86) 

* 

The  cross-terms  resulting  from  the  product  of  X (w)  and  X (w) 

P P 

are  considered  to  be  zero  since  it  was  initially  assumed  that 
X(w  + w ) and  X(w  - w ) do  not  overlap.  From  Eq  (86)  it  can 
be  seen  that  two  cases  must  be  considered  in  determining  the 
power  spectral  density,  namely,  when  n / m and  when  n = m. 
From  Eq  (73)  it  was  shown  that  E{ ( -l)-n( -l)_m]  = 0 when 
n / m.  Therefore,  substituting  this  result  in  Eq  (86) 

E[|Xp(w)|2}  = 0 (87) 

when  n / m.  Similarily,  substituting  the  result  of  Eq  (74) 
into  Eq  (86)  results  in 

- N N 9 
E(|Xfw)| 2}  = E E T2/4 
F n=  o m=  o 


sin  [(w  + w )T/2] 
[(w  + wc)T/2]2 
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+ 


(88) 


when  n 


p 

sin  [(w  - wc)T/2] 

[(w  - wc)T/2]2 

m.  However,  since  m = n,  Eq  (88)  can  be  reduced  to 


E{|X  (w)|2}  = 


2 

sin  [(w  + wc)T/2] 
[(w  - wc)T/2]2 


+ 


sin2  [(w  - wc)T/2] 
[(w  - wc)T/2]2 


(89) 


Thus,  the  power  spectral  density  of  the  output  signal  of  a 
BPSK  modem  is  given  by  Eq  (89). 

QPSK  Modem  Characterization 

Using  the  modem  model  shown  in  Figure  4 of  Chapter  III 
and  recalling  that  a QPSK  modem  can  be  thought  of  as  two 
BPSK  modems  in  parallel,  then  the  output  of  a QPSK  modem  is 
represented  as 


xp(t)  = xpl(t)  + xp2(t) 


(90) 


where 


N b 

x (t)  = Z (-1)  n cos  w tP(t  - 2nT) 
J n=o 


and  is  the  in-phase  component  of  2^p(‘t)» 

N b 

3C_ 2 ( "k ) = E (-1)  sin  w tP(t  - 2mT) 
p m=o 


(91) 


(92) 
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is  the  quadrature  component  of  xp(t).  The  factor  of  2 in 
the  terms  P(t  - 2nT)  and  P(t  - 2inT)  of  Eqs  (83)  and  (84), 
respectively,  results  from  the  forming  of  the  two  bit 
streams,  x^p(t)  and  Xp2(t),  from  the  binary  input  signal. 

The  signaling  interval  for  the  bit  streams,  x ^(t)  and  Xp2(t) 
is  double  that  of  the  input  signaling  interval  (Ref  1.3:33 6). 
Since  Xp^CiO  is  independent  of  Xp^^K  'the  mean  value,  auto- 
correlation, and  power  spectral  density  of  each  of  these 
processes  can  be  calculated  independently.  Since  Eq  (91) 
is  similar  to  Eq  (64),  with  the  only  difference  being  the 
factor  of  2 in  the  term  P(t  - 2nT)  in  Eq  (58).  then  the  mean 
value  and  autocorrelation  calculated  for  Xp(t)  in  the  section 
on  BPSK  Modem  Characterization  is  the  same  for  xpl(t)  since 
the  term  P(t  - 2nT)  is  a constant  for  these  calculations. 
Additionally,  Eq  (92)  differs  from  Eq  (66)  only  in  the  fac- 
tor of  2 in  the  term  P(t  - 2mT)  and  in  the  carrier  term  of 

sin  w t.  However,  in  the  calculation  of  the  mean  value,  the 

c 

carrier  term  is  constant,  therefore,  the  mean  value  for  xp2 
(t)  is  that  of  Eq  (70).  Therefore,  the  mean  of  xp(t)  in 
Eq  (90)  is  zero  since  the  mean  of  x^^t)  and  Xp^'t)  are  ^oth 
zero . 

The  autocorrelation  function  of  Eq  (92)  is  the  same  as 
that  of  Eq  (?2)  with  the  term  sin  wQt1  sin  wct2  substituted 
in  for  cos  wct^  cos  wct2>  Therefore,  the  autocorrelation 
function  of  Eq  (92)  is 
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f 


0 when  tp  and  t2  are  not  in 

the  same  signaling  interval 

sin  w t,  sin  w t0  when  tn  and  t0  are  in  the 
cl  c c 1 c. 

same  signaling  interval 

(93) 


Therefore,  the  autocorrelation  function  of  xp(t)  in  Eq  (90) 
is 


H( tx,t2) 


0 when  tp  and  t2  are  not  in 

the  same  signaling  interval 

cos  w t^  cos  wct2  w^en  and  are  in  "the 

.+  sin  wc‘tj_  sin  wct2  same  signaling  interval 

(94) 


which  indicates  that  Xp(t)  is  a nonstationary  process  in  a 
wide-sense.  As  a result,  Eq  (44)  is  used  to  find  the  power 
spectral  density  which  is  represented  as 

E[|Xp(w)|2}  = E{Xp(w)X*(w)} 

= E{[Xpl(w)  + Xp2(w) ] [Xpl(w)  + Xp2(w)]}  (95) 

Eq  (95)  is  found  by  transforming  Eq  (90)  from  the  time- 
domain  to  the  frequency-domain  and  multiplying  this  result 
by  its  complex  conjugate.  Eq  (95)  is  reduced  to 


E[|Xp(w)|2}  = E{Xpl(w)X*1(w)  + Xp2(w)X*2(w)}  (96) 

and  the  product  of  cross  terms,  Xp2(w)Xpl(w)  and  x;2<w)xpl(w) 
are  considered  zero  as  a result  of  assuming  that  the  negative 
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frequency  terms  and  positive  frequency  terms  do  not  overlap. 
Using  the  definition  of  expectation  of  a sum  of  two  inde- 
pendent processes,  Eq  (93)  is  represented  as 

E{|Xp(w)|2}  = E{Xpl(w)X*1(w)}  + E{Xp2(w)X*2(w)}  (97) 

The  power  spectral  density  of  the  first  term  in  Eq  (97)  is 
similar  to  that  calculated  in  the  section  on  BPSK  Modem  Char- 
acterization. This  can  be  seen,  since  the  transformation  of 
Eq  (91)  is  that  of  Eq  (82)  with  T = 2T . Making  the  neces- 
sary substitutions  and  using  Eqs  (84),  (85)»  and  (86),  the 
power  spectral  density  of  the  first  term  in  Eq  (97)  is  equal 
to  that  of  Eq  (89)  with  T = 2T . 

The  frequency-domain  transformation  of  Eq  (92)  is  found 
by  using 


F [sin  w t] 


1 

2 


6 (w  + wc) 


(w  - wQ) 


(98) 


and  multiplying  this  result  with  that  of  Eq  (81)  with  T = 2T 
N 

and  the  term  £ (-1)  . The  transformation  of  Eq  (92)  is 

n=o 


then 


(w) 


N b 
E (-D  m(  T ) j 
n=o 


- jT( w+w  ) ( 2n+l) 

sin  [(w  + w )T]  e 
(w  + wc)T 


-jT(w-w  )(2m+l)  ' 

sin  [ (w  - wc)T]  e 
(w  - wc)T 


(99) 


Using  the  procedure  outlined  in  the  section  on  BPSK  Modem 
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Characterization  on  finding  the  power  spectral  density,  the 
second  term  in  Eq  (97)  has  power  spectral  density  equal  to 
that  of  the  first  term.  Therefore,  the  power  spectral  den- 
sity of  Xp(t)  in  Eq  (90)  is 


2 

= NT 


+ NT2 


2 

- 2NT 


sin2  [ (w 

+ wq)T] 

sin2  [(w  - wc)T] 

i — i 

+ 

wc)T]2 

[(w  - wc)T]2 

sin2  [ (w 

+ wq)T] 

sin2  [(w  - wc)T] 

[(w  + wc)T]2 

C(w  - wc)T]2 

sin2  [(w 

+ WC)T] 

p 

sin  [(w  - w )T] 

[(w  + wc)t]2  [(w  - wc)T]2  J 


(100) 


BFSK  Modem  Character izat ion 

Using  the  modem  model  shown  in  Figure  6 of  Chapter  III 
with  an  input  a random  pulse  train,  the  output  signal  of  a 
BFSK  modem  is  represented  as 


N b 

x ( t)  = 2 cos  (w  + (-1)  nAw)tP(t  - nT)  (101) 

p n=  o 


where 


b 

— n 


random  variable  equally  likely  to  be  one 
or  zero 


AW 


frequency  deviation  from  wc 


N 

E P(t  - nT) 
n=o 


the  n^  signaling  interval  of  the  input 
signal  summed  over  n intervals  with  each 
pulse  having  period  T 


w 


c 


carrier  frequency 


48 


The  mean  or  expected  value  of  xp(t)  is 


E[xp(t)}  = E 


N k, 

E P(t  - nT)  cos  (w  + (-1)  Aw)t 
n=o 


N 

E P(t  - nT)E 
n=o 


cos  (wc  + ( -l)  nAw) t 


(102) 


Using  a trigonometric  identity  to  expand  the  cosine  argu- 
ment, (i.e.,  cos  (u  + v)  = cos  u cos  v - sin  u sin  v ) 
Eq  (102)  reduces  to 


Ety 


N b 

x^t)}  = E P(t  - nT)E{cos  w t cos  (-1)  nAwt 
n=c  c 


— n. 


- sin  w t sin  (-1)  Awt} 


N 


= E P(t  - nT)  [(cos  w t)E{ cos  (-1)  ntAw} 
n=  o 1 


b 

- (sin  wct)E{sin  (-1)  nAwt}  (103) 


Now  substituting  in  the  values  of  zero  and  one  for  bn  and 
evaluating,  Eq  (103)  reduces  to 


N 

x^t)}  = E P(t  - nT)  COS  W t COS  Awt 
n=o  c 


E[Xp< 


(104) 


which  is  the  mean  value  of  x (t) 

~P 


Now  the  autocorrelation  function  of  Xp(t)  using  Eq 
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( 21)  is 


R(ti,t2)  = E{xp(t1)xp(t2)} 

= E{E  E P(t1  - nT)P(t2  - mT)  cos  (wc 

b b 

+ (-1)  nAw)t1  cos  (wc  + (-1)  mAw)t2} 

N N 

= E E PCt-j-  - nT)P(t2  - mT)E{cos  (wq 
n=  o m=  o 

b b 

+ (-1)  nAw)t1  cos  (wc  + (-1)  mAw)t2}  (105) 

In  Eq  (105),  xp(t)  is  sampled  in  the  same  manner  as  is  xp(t) 
in  Eq  (72)  . 

The  BFSK  autocorrelation  function,  Eq  (105),  is  deter- 
mined using  a procedure  similar  to  that  used  in  determining 
the  BPSK  autocorrelation  function,  Eq  (72).  In  determining 
the  autocorrelation  function  of  Eq  (105),  two  cases  are  con- 
sidered. The  first  case  is  when  the  sample  times,  tp  and  t2# 
are  not  in  the  same  signaling  interval,  n / m.  The  second 
case  is  when  the  sample  times,  tp  and  t2,  are  in  the  same 
signaling  interval,  n = m.  With  the  use  of  Tables  I and 
II,  the  definition  and  properties  of  expectation  and  the 
procedures  for  evaluating  Eq  (73) » the  autocorrelation  func- 
tion of  Eq  (105)  is  given  as 

R(t1,t2)  = 0 (106) 

when  n / m and 
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(107) 


R(tltt2)  = sin  wct-j_  sin  wct2  s^n  Aw_ti  sin  Aw_t2 


+ cos  w t,  cos  w t0  cos  Awt,  cos  AWt0 
c i c 2 1 2 


when  n = m.  Combining  Eqs  (106)  and  (107) 


cos  w t,  cos  w t0  cos  Awt,  when  t,  and  t0  are 
cl  c 2 1 1 2 


cos  Awt^  in  the  same  sig- 


+ sin  w t^  sin  wc‘t2  s^n  Aw_tp  naling  interval 


R (tlft2) 


sin  Awt , 


when  tp  and  t2  are 
not  in  the  same 
signaling  interval 
(108) 


Since  Eq  (108)  indicates  that  Xp(t)  is  not  a wide  sense 
stationary  process,  Eq  (44)  is  used  to  determine  the  power 
spectral  density.  The  power  spectral  density  of  xp(t)  is 
represented  as 


E{|Xp(w)l  } = E{Xp(w)X*(w)} 


(109) 


where  X (w)  is  the  Fourier  transform  of  Eq  (101)  in  the  fre- 
quency-domain. X*(w)  is  the  complex  conjugate  of  X (w) . 

P bn  P 

With  wn  = (w  + (-1)  aw)  in  Eq  (101),  the  Fourier  trans- 


form of  Xp(t)  is  represented  as 


X (w)  = z T/2 


- j T ( w+w  ) ( n+ 1/2 ) 

sin  [(w  + w )T/2]  e 


> + wn)T/2 


V 


+ 


(110) 


sin  [(w  - w )T/2]  e 
(w  - wn)T/2 


-jT(w-w) (n+l/2) 


-n 


Substituting  Eq  (110)  and  dts  complex  conjugate  into  Eq 
(109),  the  power  spectral  density  of  Xp('t)  Is 


, N N 9 
E{|X(w)l 2]  = £ £ T2/4 

y n=o  m=o 


sin  [(w  + wn)T/2] 
(w  + wn)T/2 


- j T ( w+ w ) ( n+ 1/2 ) - j T ( w+ w ) ( m+ 1/2 ) 
sin  [(w  + wm)T/2]  e e 

(w  + wm)T/2 


sin  [(w  - w )T/2] 
(w  - wn)T/2 


-jT(w-w) (n+l/2)  - jT(w-w  ) (m+l/2) 

sin  [(w  - wn)T/2]  e e 


(w  - wm)T/2 


(111) 


The  cross-terms  resulting  from  Eq  (109)  are  assumed  zero 
and  are  based  on  a similar  assumption  made  in  deriving  Eq 
(86),  Again,  two  cases  must  be  considered  in  determining 
the  power  spectral  density,  namely,  when  n / m and  n = m. 
After  performing  the  necessary  mathematical  manipulations, 
the  power  spectral  density  of  Xp("t)  Is 


T^N 

E[|X(w)|2}  = “ “ 


sin  [(w  ± WpjT/2] 
[(w  ± wx)t/2]2 
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sin2  [(w  ± w2)T/2] 
[(w  ± w2)t/2]2 


(112) 


where 


= wc  + aw 

w2  = wc  - AW 

The  carrier  frequencies,  w-^  and  w2,  are  those  indicated  in 
Eqs  (52)  and  (53)*  respectively. 

DBFSK  Modem  Characterization 

Using  the  modem  model  shown  in  Figure  9 of  Chapter  III 
with  the  input  a random  pulse  train  encoded  as  a duobinary 
signal,  the  output  of  a DBFSK  modem  is  represented  as 


x_(t)  = E h(t)P(t  - nT)  cos  (wc  + xdb(t)Aw)t  (113) 

*•  n — rs  ^ 


where 


h(t) 


= the  duobinary  transformation 

= random  process  having  value  of  1,  0, 
and  -1,  but  not  equally  likely 


n ^ 

E P(t  - nT)  = the  n signaling  interval 


w = carrier  frequency 

From  information  published  on  duobinary  encoding  (Refs 
8; 214-218,  11:126-130)  and  digital  transmission  (Ref  2:1501- 
1511),  the  process  Xp(t)  is  nonstationary.  Therefore,  the 
power  spectral  density  is  found  using  Eq  (44) . Using  Eq 
(81),  the  Fourier  transform  of  P(t  - nT);  Eq  (54),  the 
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Fourier  transform  of  h(t);  and  the  Fourier  Transform  of 
cos  [w  + xdb(  t)Aw]t  ; the  Fourier  transform  of  Eq  (113) 
is 


Xp(w) 


N 

E Vi/!  (T/2) 
n=o 


cos  [(w  + w )T/2] 


- j T ( w+w  ) ( n+ 1/2 ) 

sin  [(w  + wn)T/2]  e 
(w  + wn)T/2 

-jT(w-w  ) (n+l/2) " 

cos  [(w  - wn)T/2]  sin  [ ( w - wn)T/2]  e 
+ (w  - wn)T/2 

(114) 

where  the  substitution,  wn  = wQ  + xdb(t)Aw  , in  Eq  (113) 
was  made  prior  to  taking  the  Fourier  transform  of  cos 
£w  + xdb(t)Aw]t  . Substituting  into  Eqs  (83)  and  (84), 
the  power  spectral  density  is 


EC|Xp(w)|2}  = E 


N N 9 

E E T78 

n=o  m=o 


cos  [(w  + wn)T/2] 


cos  [(w  + wm)T/2]  sin  [(w  + wn)T/2] 

(w  + wn)T/2 

-jT(w+w  ) (n+l/2)  -jT(w+w  ) (m+l/2) 
sin  [(w  + wm)T/2]  e e 

(w  + wm)T/2 
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cos  [(w  - wn)T/2]  cos  [(w  - wm)T/2]  sin  [(w  - wn)T/2] 


(w  - w^T/2 


sin  [(w  - wm)T/2]  e 
(w  - wm)T/2 


- jT(w-w  ) (n+1/2)  -jT( w-w  ) (m+l/2) 


(115) 


Again,  the  cross-terms  are  assumed  zero.  With  wn  = wc 

+ xdb(t)Aw  and  with  xdb(t)  having  values  of  +1,  0,  and  -1, 

then  w = w + aw  when  xH-h(t)  ~ +1  ; w = w.  when 

c ci  d n c 

xdt>(  = 0 ; 3X1(1  wn  = wc  " Aw  when  xdb^^  = -1  ' 

Therefore,  the  three  values  of  xdb(t)  will  result  in  w^ 

having  three  values.  The  probability,  P,  of  the  occurrence 

of  the  levels  of  the  duobinary  signal  are  given  as  P(+l)  = 

P(-l)  = 1/4  and  the  P(0)  = 1/2  (Ref  8:218). 

From  Eq  (115)  it  can  be  seen  that  two  cases  must  be 

considered  in  determining  the  power  spectral  density,  namely, 

when  n / m and  when  n = m.  After  performing  the  necessary 

mathematical  manipulations,  the  power  spectral  density  of 


Xp(t)  i: 


2 nt2  sin2  [(w  ± w )T/2] 

E{|X  (w)|2}  = ^4-  — 5— 

P [(w  ± w-^T/2]2 


2 sin2  [(w  ± w )T/2] 
[(w  ± w?)T/2]2 
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+ 


sin2  [(w  ± w2)t/2] 
[(w  ± w2)T/2]2 


(116) 


where 

w-^  = w£  + aw 
W2  = Wc 

W0  = W„  - AW 

3 c 

The  carrier  frequencies,  w^,  w2,  and  w^,  are  those  indicated 
in  Eqs  (57),  (58),  and  (59)- 


BASK  Modem  Characterizat ion 

Using  the  modem  model  shown  in  Figure  10  in  Chapter  III 
with  an  input  a random  pulse  train,  the  output  signal  of  a 
BASK  modem  is  represented  as 


xp(t) 


N X* 

1 + Z (-1)  nP(t  - nT) 
n=o 


cos  wct 


(117) 


The  mean  or  expected  value  of  x(l0  is 


N b 1 

CXp(t)}  = E 

1 + Z (-1)  ^>(t  - nT) 
n=o 

cos  w t 
c 

N 6 

= cos  w t + Z P(  t - nT)  cos  w t Ef(-l)  j (118) 
n=o 

bn 

From  Eq  (69),  E{(-1)  } = 0 , therefore,  the  expected  or 

mean  value  of  Xpt^)  is 

(119) 


E{25p(‘t)}  = cos  wct 
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The  autocorrelation  function  of  xp(t)  using  Eq  (21)  is 


R(tltt2)  = E{xp(t1)2cp(t2)} 


N V 

E [1  + E (-1)  T(t-.  - nT)~l  cos  w t, 

n=o  x c x 


N b 


[1  + E (-1)  mP(t2  - mT)]  cos  wct2 


N N 


= cos  WQt1  cos  wct2  + E E P(t1  - nT)P(t2  - mT) 


n=o  m=o 


b b 

cos  w t,  cos  w t„  E{(-1)  n(-l)  m} 


(120) 


Using  the  results  of  Eqs  (73)  and  (74),  the  autocorrelation 


function,  Eq  (120),  is 


J 


+ 


f 


I 


t 


N b 
£ (-1)  m(T/2) 
n=o 


% . -jT(w+w  )(n+l/2) 

sin  [(w  + w )T/2"]  e 

(w  + wc)T/2 


N b 

+ £ (-1)  n(T/2) 

n=o 


-jT(w-w  ) (n+l/2) 

sin  [(w  - wc)T/2]  e 
(w  - wq)t/2 


(122) 


where  the  term  1/2  [s(w  + w ) + 6(w  - w )]  is  the  Fourier 
transform  of  product  of  the  carrier  signal,  cos  wpt,  and  the 
dc  bias,  1. 

Substituting  Eq  (122)  and  its  complex  conjugate  into 
Eq  (44)  and  using  the  results  from  Eqs  (86),  (87),  and  (88), 
the  power  spectral  density  is  given  as 

E[|X_.(W)|2}  = 1/4  6 (w  + w ) + 1/4  5 (w  - w ) 


sin^  [(w  + w )T/2] 


[( w + wq)T/2 


2 


sin2  [(w  - wc)T/2] 
[(w  - wc)T/2]2 


(123) 


The  firsv-  two  terms  in  Eq  (123)  are  delta  functions  located 

at  the  carrier  frequencies,  w and  -w  , and  have  weight  1/2. 

c c 

DSBAM  Modem  Characterization 

Using  the  modem  model  shown  in  Figure  11  of  Chapter  III, 
and  with  an  input  signal  given  by  Eq  (65),  the  output  signal 
of  a DSBAM  modem  is  represented  as 

x(t)  = [1  + a(t)]  cos  wct 

= (y  cos  wst  + 1)  cos  w£t  (124) 
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where  a(t)  = y cos  wct  is  the  input  signal.  The  Fourier 

s 

transform  of  Eq  (124)  is  given  as  (Ref  12:423) 

X(w)  = J 6 (w  - ws  + wc)  + J 6 (w  + wg  + wc) 

+ I 6 (w  - wc  - ws)  + l 6 (w  + wg  - wc) 

+ | (w  ~ wc)  + i (w  + Wc)  (125^ 

Using  Eqs  (34)  and  (125),  the  power  spectral  density  or  more 
correctly  the  power  spectrum  is  given  as 

2 

E[|Xp(w)|2}  = | (w  - wc)  + | (w  + wc)  + 6 (w  - ws  + wc) 

2 2 
+ 16  6 (w  + ws  + wc)  + 16  6 (w  - wc  - ws) 

Y2 

+ 6 (w  + wg  - wc)  (126) 

where  the  terms  6/4(w  - w ) and  6/4(w  + w ) are  delta  func- 

C 4/ 

tions  located  at  the  carrier  frequencies,  w and  -w  , with 
each  having  weight  1/4.  The  remaining  terms  are  delta  func- 

p 

tions  having  weight  y /l6,  where  0 < y < 1. 

DSBAM-SC  Modem  Characterization 

Using  the  modem  model  shown  in  Figure  12  of  Chapter  III, 
and  with  an  input  signal  given  by  Eq  (65),  the  output  signal 
of  a DSBAM-SC  modem  is  represented  as 

x(t)  = a(t)  cos  w t 


59 


= Y cos  W3  COS  wfit 


(127) 


where  a(t)  = y cos  w~  is  ‘the  lnPu_t  signal. 

o 

The  Fourier  transform  of  Eq  (127)  is  given  as  (Ref  12: 

428) 

Xp(w)  = l 6 (w  - ws  + wc)  + * 6 (w  + ws  + wc) 

+ 4 6 (w  - wg  - wc)  + ^ 6 (w  + wg  - wc)  (128) 

Since  Eq  (128)  differs  from  Eq  (125)  in  that  the  carrier 
term,  [6/2(w  - w ) + 6/2(w  + w )]  has  been  suppressed, 
the  power  spectral  density  or  more  correctly  the  power 
spectrum  is  given  as 

2 

E[|Xp(w)|2}  = [6  (w  - wc  - wg)  + 6 (w  + wc  + ws) 

+6  (w  - wc  + wg)  +6  ( w + wc  - ws)]  (129) 
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V . Results  and  Discussion 


This  chapter  is  divided  into  two  sections  . The  first 
section  covers  the  characterization  results  and  the  second 
covers  the  data  rate  results. 

Modem  Identification  Results 

The  primary  characterization  of  all  modems  presented 
in  this  study  is  the  power  spectral  density  of  the  output 
signal  from  each  modem.  Table  III  lists  the  modems  char- 
acterized and  the  equation  for  the  power  spectral  density 
of  the  output  signal  from  each  modem.  Comparing  the  equa- 
tions in  Table  III,  it  can  be  seen  that  five  of  the  seven 
modems  listed  can  be  individually  identified  by  the  power 
spectral  density.  Only  the  PSK  modems  (BPSK  and  QPSK) 
cannot  be  distinguished  between  themselves  as  the  only 
difference  in  their  power  spectral  density  equations  is  a 
constant  times  the  sine  squared  terms  and  the  width  of 
the  signaling  interval,  T.  Neither  of  these  differences 
can  be  used  to  distir guish  the  BPSK  modem  from  the  QPSK 
modem.  Thus,  Figure  13  represents  a decision  tree  scheme 
for  identifying  the  modems  listed  in  Table  III  using  as 
the  decision  criteria,  the  power  spectral  density. 

Since  the  BPSK  and  QPSK  modems  are  distinguished  only 
in  the  general  category  of  PSK  modems,  a second  character- 
ization was  chosen  to  distinguish  between  individual  modems 
in  this  general  category.  The  characterization  chosen  to 
differentiate  the  BPSK  modem  from  the  QPSK  modem  was  the 
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TABLE  III 

Characterization  of  Modems  by  Power  Spectral  Density 


Modem 


Power  Spectral  Density,  E{|XWI2}  (dBw) 


BPSK 


,2  sin2  [(w  + w )t/2]  sin2  [(w  - w JT/2] 
— + 

4 [(w  + w )T/2]2  [(w  - w )T/2]2 

. c c 


QPSK 


„ [ sin2  [(w  + w )T]  sin2  [(w  - w )T] 
2NT  5 + — 


[(w  + wc)T]2 


[(w  - wc)T]2 


BFSK 


nt2  [ sin2  [(w  ± w1)T/2]  sin2  [(w  ± wg)T/2] 


C(w  ± w1)T/2]: 


[(w  ± w2)T/2]: 


Tyjm 2 sin2  [(w  ± w,  )t]  2 sin2  L(w  ± w?)T] 

DBFSK  £4-  J 

[(w  ± W-^T]2  [(w  ± w2)T]2 


sin2  [ ( w ± w^)T] 
[(w  ± w^)T]2 


n„2  sin  [(w  ± w )T/2]  , 

BASK  5 — — + 4 [6(w  + w ) 

4 [(w  ± wc)T/2]2  4 

+ 6 ( w - wc)] 


DSBAM 


1 

- [s(w  + wc)  + 6 ( w - wc)]  + [6(w  + wc  + wg) 


+ 6 (w  - wq  - wg)  + 5(w  - wc  + wg) 


+ 6 (w  + wc  - ws)] 


DSBAM-SC  [6(w  - wq  - wg)  + 6(w  + wg  + wg) 


+ 6 ( w - wc  + ws)  + 6 ( w + wc  - ws)] 


Figure  13 . Modem  Identification  Decision  Tree  Scheme 
Using  Power  Spectral  Density  as  the 
Decision  Criteria 


autocorrelation  function.  While  the  power  spectrum  and  the 
correlation  form  a Fourier  transform  pair  for  a wide-sense 
stationary  process,  it  is  not  so  for  nonstationary  processes, 
and  the  processes  for  BPSK  and  QPSK  are  nonstationary. 

The  autocorrelation  functions  calculated  in  Chapter  IV 
for  the  BPSK  and  QPSK  modems  are  listed  in  Table  IV.  Com- 
parison of  these  functions  shows  that  the  autocorrelation 
function  can  be  used  to  distinguish  between  a BPSK  modem  and 
a QPSK  modem.  Figure  14  represents  a decision  tree  scheme 
for  identifying  all  the  modems  listed  in  Table  III.  The 
term  represents  the  power  spectral  densities  of  all  the 
modems  being  analyzed.  In  area  HQ  of  Figure  14,  the  cri- 
teria used  to  make  the  decision  in  identifying  a particular 
modem  is  the  power  spectral  density.  In  the  area  represented 
by  H-^,  the  decision  is  based  on  the  autocorrelation  function. 
However,  for  the  autocorrelation  function  to  be  a valid  means 
of  identifying  the  type  of  modem,  the  sampling  of  the  output 
signal,  Xp("t)»  must  occur  during  the  same  signaling  interval 
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TABLE  IV 

Characterization  of  Modems  by  Autocorrelation  Function 


Modem 


Autocorrelation  Function,  R(t^,t2 


0 when  t^  and  t2  are  not  in 

the  same  signaling  interval 

cos  w t,  cos  w t„  when  t,  and  t„  are  in  the 
cl  c 2 1 2 

same  signaling  interval 

0 when  t-^  and  t2  are  not  in 

tie  same  signaling  interval 

sin  w t,  sin  w t0  when  t,  and  t„  are  in  the 
cl  c 2 1 2 

+ cos  w t^  cos  wct2  same  signaling  interval 


Figure  14.  Modem  Identification  Decision  Tree  Scheme 
Using  Power  Spectral  Density  and  the  Auto- 
correlation Function  as  the  Decision 
Criteria 
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A correlator  and  frequency  spectrum 


for  times  t-^  and  t^. 
analyzer  can  be  used  to  determine  the  autocorrelation  func- 
tion and  power  spectral  density,  respectively. 

Data  Rate  Identification  Results 

The  equations  for  the  power  spectral  densities  of  the 
modems  listed  in  Table  III  (except  for  the  DSBAM  and  DSBAM- 
SC  modems),  were  used  to  construct  power  spectral  density 
curves  for  various  data  rates.  In  the  equations  listed  in 
Table  III  (except  for  the  QPSK  modem),  the  data  rate  is  rep- 
resented as  T~\  where  T is  the  signaling  interval.  In  the 
equation  for  the  QPSK  modem,  T_1  represents  symbol  rate  and 
T is  the  signaling  interval.  The  pcwer  spectral  curves 
(only  the  positive  half)  generated  are  presented  in  Appen- 
dix A,  Figures  15  through  Jl,  along  with  a list  of  the  data 
rates  and  carrier  frequencies  used  in  constructing  the 
curves.  The  data  rates,  given  in  bits  per  second  (BPS), 
and  the  carrier  frequencies,  given  in  hertz  (Hz),  are  valid 
modem  specifications. 

The  power  spectral  density  curves  in  Appendix  A were 
used  to  determine  the  data  rate  (symbol  rate  for  the  QPSK 
modem)  of  the  output  signal,  Xp(t),  primarily  by  using  two 
different  methods. 

The  first  method  involved  measuring  the  width  of  the 
main  lobe  (lobe  with  the  largest  amplitude^  at  the  crossover 
points  on  the  frequency  axis  and  dividing  this  value  by  two. 
For  the  QPSK  modem,  this  method  gives  the  symbol  rate.  For 
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the  BFSK  modem,  Figures  23,  24,  and  25,  only  one  of  the 
main  lobes  was  used,  and  for  Figure  26,  the  entire  curve 
was  assumed  to  be  the  main  lobe. 

The  second  method  involved  measuring  the  width  of  the 
main  lobe  at  the  3 dB  points.  For  the  BFSK  modem,  Figures 
23,  24,  and  25,  only  one  of  the  main  lobes  was  used,  and 
for  Figure  26,  the  entire  curve  was  assumed  to  be  the  main 
lobe . 

The  analytical  results  of  these  two  methods  are  listed 
in  Table  V.  Using  the  main  lobe  crossover  method  proved  to 
be  a valid  method  for  determining  the  data  rate  except  for 
the  BFSK  modem  at  1200  BPS  and  the  DBFSK  modem  at  2400  BPS. 
The  analytical  results  of  the  3 dB  point  method,  Table  V, 
were  not  as  valid  as  the  results  for  the  main  lobe  cross- 
over method,  especially  for  the  BFSK  modem  at  1200  BPS. 

A secondary  or  third  method  was  tried  to  more  accu- 
rately determine  the  data  rate  of  the  BFSK  modem  at  1200 
BPS.  The  assumption  used  in  the  first  two  methods,  that 
the  power  spectral  density  curve  for  the  FBSK  modem  at  1200 
BPS,  Figure  26,  was  one  single  main  lobe,  was  changed.  The 
curve  in  Figure  26  was  now  assumed  to  be  composed  of  two 
main  lobes  which  if  Figures  23,  24,  and  25  are  viewed,  sup- 
ports this  assumption.  As  the  data  rate  is  increased,  the 
main  lobes  get  closer  together  until  they  overlap  or  merge 
as  is  shown  in  Figure  26  when  the  data  rate  is  1200  BPS. 
Thus,  the  third  method  was  based  on  measuring  the  width 
between  the  two  peaks  of  curve  shown  in  Figure  26.  As  a 
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TABLE  V 

Predicted  Versus  Actual  Modem  Data  Rate  Values 


Modem 


BPSK 


QPSK 


BFSK 


DBFSK 

BASK 


Data  Rate  (BPS)* 

Actual 

3 dB  Point  Method 

Main  Lobe 
Crossover 
Method 

15 
30 
60 
1200 

150 

30Q 

600 

1200 

150 

300 

600 

1200 

2400 

150 

300 

600 

1200 


120 

270 

540 

1060 

120 

280 

580 

2280 

1530 

120 

270 

540 

1060 


1200 

150 

300 

600 

1200 

150 

300 

600 

1470 

1500 

150 

300 

600 

1200 


For  the  QPSK  modem,  the  rates  listed  are  symbol  rates 
expressed  as  symbols  per  second  (SPS). 


result,  this  method  produced  a data  rate  of  1200  BPS  as 
compared  to  the  actual  data  rate  of  1200  BPS. 

A fourth  method  was  tried  to  more  accurately  determine 
the  data  rate  of  the  DBFSK  modem  at  2400  BPS.  From  the  dis- 
cussion on  the  DBFSK  modem  in  Chapter  III  and  from  (Ref  11: 
128),  the  major  advantage  of  duobinary  encoding  is  the 
pseudo-compression  of  the  spectrum  along  the  frequency  axis 
when  the  input  binary  signal  is  in  the  form  of  a rectangular 
pulse  train.  This  means  that  certain  parameters  in  the 
spectrum,  which  are  sometimes  used  to  define  bandwidth  in 
a loose  sense,  are  also  halved.  One  of  these  is  the  fre- 
quency at  which  the  spectrum  goes  to  zero  or  the  first 
crossover  point  of  the  spectrum  on  the  frequency  axis.  In 
essence,  duobinary  encoding  allows  for  higher  data  rates  to 
be  transmitted  in  the  same  bandwidth  (in  loose  terms)  as 
lower  data  rates  of  binary  signals.  The  binary  signal  is 
restricted  in  this  case  to  a rectangular  pulse  train. 
Therefore,  the  original  method  used  to  obtain  the  result  in 
the  main  lobe  crossover  method  (Table  V)  was  changed  in  that 
only  the  width  of  the  main  lobe  at  the  crossover  points  on 
the  frequency  axis  is  used  as  the  measure  of  the  data  rate . 
As  a result,  this  method  produced  a data  rate  of  3000  PBS 
as  compared  to  the  actual  data  rate  of  2400  BPS. 

Of  the  methods  used  to  determine  the  data  rate  (symbol 
rate  for  QPSK  modem) , only  the  main  lobe  crossover  method 
was  successful  except  for  determining  the  data  rates  for 
the  BFSK  modem  at  1200  BPS  and  the  DBFSK  modem  at  2400  BPS. 
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This  method  is  not  a valid  predictor  of  the  data  rate  for 
these  modems  at  the  data  rate  specified  because  of  the  un- 
certainty of  the  width  of  the  main  lobes  due  to  the  merger 
of  the  two  lobes  in  the  power  spectral  curve  for  the  BFSK 
modem  and  of  the  merger  of  the  three  lobes  in  the  powar 
spectral  curve  for  the  DBFSK  modem.  An  alternate  method 
was  used  to  determine  the  data  rate  of  the  BFSK  modem  when 
the  two  main  lobes  merged  and  although  it  proved  success- 
ful, it  is  considered  valid  only  in  the  case  used.  The 
alternate  method  used  to  determine  the  data  rate  of  the 
DBFSK  modem  proved  unsuccessful. 


69 


VI.  Conclusions  and  Recommendations 


Conclusions 

The  major  conclusion  reached  in  this  study  is  that 
modems  can  be  characterized  by  the  real  time  measurement  of 
the  power  spectral  density  and  the  autocorrelation  function 
of  the  output  signal  from  the  modem.  This  characterization 
in  turn  leads  to  the  identification  of  the  modulation  scheme 
being  used  as  each  modem  or  group  of  modems  has  a unique 
power  spectral  density.  To  identify  individual  modems  with- 
in a group,  the  autocorrelation  function  can  be  used. 

Another  conclusion  is  that,  in  most  instances,  the  data 
rate  (symbol  rate  for  QPSK  modem)  can  be  determined  by  ob- 
serving the  width  of  the  main  lobes  of  the  power  spectral 
curves.  However,  this  is  valid  only  when  the  main  lobes  do 
not  merge  as  in  the  BFSK  and  DBFSK  modems. 

Recommendations 

The  power  spectral  density  curves  presented  in  Appendix 
A could  be  used  as  templates  or  as  overlays  to  be  placed  on 
spectrum  analyzers  in  order  that  an  operator  might  determine 
the  modulation  scheme  of  the  modem.  Additionally,  the  ope- 
rator could  determine  the  data  rate  (symbol  rate  for  QPSK 
modem)  at  which  the  modem  is  operating  by  measuring  half 
the  width  of  the  main  lobe  at  the  crossover  points  on  the 
frequency  axis  of  the  power  spectral  density  curve.  How- 
ever, this  is  a valid  test  only  when  main  lobes  do  not  merge 
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as  in  the  case  for  the  BPSK  modem  at  1200  BPS  and  the  DBFSK 
modem  at  2400  BPS. 

As  a result  of  the  conclusions  and  the  research  re- 
quired for  this  study,  the  following  recommendations  are 
submitted  for  future  studies  in  this  area. 

1.  Since  this  study  presents  only  theoretical  results 
on  the  characterization  of  modems  by  the  power 
spectral  density  and  the  autocorrelation  func- 
tions, empirical  data  is  needed  to  test  the  con- 
clusions. Additionally,  the  study  should  include 
how  additive  interference  (i.e.,  additive  white 
gaussian  noise)  and  other  multiplicative  effects 
such  as  fading  (a  time-varying  channel  attenuation 
and  phase  shifting  of  the  desired  wave)  affect 
the  characterization  of  modems  using  second  moment 
characterizations . 

2.  Additional  research  is  needed  to  find  a method 
for  determining  the  data  rate  at  which  a modem 
is  operating. 

3.  A future  study  would  be  to  investigate  the  possi- 
bility of  using  some  feature  or  parameter,  once  a 
modem  has  been  identified,  to  do  degradation  anal- 
ysis on  the  received  signal  in  a technical  control 
facility.  Developing  such  a method  could  possibly 
aid  in  assessing  the  performance  of  communications 
links  and  enable  action  to  be  taken  before  serious 
degradation  occurs. 
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4.  Research  should  be  conducted  in  automating  the 
identification  of  modem  type  and  data  rate. 

One  approach  in  automating  the  identification  of  modem 
type  would  be  to  use  a Fast  Fourier  Transform  (FFT)  algo- 
rithm to  transform  the  modem's  output  signal  (or  received 
signal  from  the  channel)  from  the  time-domain  to  the  fre- 
quency-domain and  calculate  the  power  spectral  density.  An 
algorithm  would  then  be  used  to  compare  the  calculated  power 
spectral  density  with  other  known  power  spectral  densities 
as  a means  of  identifying  the  type  of  modem.  It  should  be 
noted  that  the  power  spectrum  and  the  correlation  form  a 
Fourier  transform-  pair  for  a wide-sense  stationary  process, 
but  not  for  nonstationary  processes.  Therefore,  taking  the 
inverse  FFT  of  the  power  spectral  density  will  not  result 
in  finding  the  complete  autocorrelation  function  if  the 
process  is  nonstationary. 
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Appendix  A 

Power  Spectral  Density  Curves 

This  appendix  contains  the  power  spectral  density 
curves  (only  the  positive  half)  for  the  modems  listed  in 
Table  III.  The  power  spectral  density  equations  (positive 
half  only)  in  Table  III  were  used  to  construct  the  curves 
for  different  values  of  T,  the  signaling  interval.  The 
inverse  of  T is  called  the  data  rate  (except  for  the  QPSK 
modem,  it  is  called  the  symbol  rate)  and  the  values  of 
T-1  and  the  carrier  frequencies  used  in  the  equations  of 
Table  III  are  listed  in  Table  VI  of  this  appendix.  The 
data  rates  are  in  bits  per  second  (BPS)  (except  for  the 
QPSK  modem,  the  rate  is  in  symbols  per  second  (SPS))  and 
the  carrier  frequencies  are  given  in  hertz  (Hz). 
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TABLE  VI 


Carrier  Frequencies  and  Data  Rates  Used  to 
Construct  Power  Spectral  Density  Curves 


Modem 

Carrier  Frequency  (Hz) 

Data  Rate  (BPS)* 

BPSK 

w = 1800 
c 

150,300,600,1200 

QPSK* 

w = 1800 
c 

150,300,600,1200 

BFSK 

wx  = 1200 
w 2 = 2400 

150,300,600,1200 

DBF  SK 

w,  = 1200 

w2  = 1800 

= 2400 

2400 

BASK 

w = 1800 
c 

150,300,600,1200 

* Rate  for  the  QPSK  modem  is  a symbol  rate  instead  of 
a data  rate . 
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Figure  17 • Power  Spectral  Density  Curve  for  BPSK  Modem  With  Data  Rate,  600  BPS 
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Figure  19*  Power  Spectral  Density  Curve  for  QPSK  Modem  With  Symbol  Rate,  150  SPS 


81 


Figure  20.  Power  Spectral  Density  Curve  for  QPSK  Modem  With  Symbol  Rate,  300  SPS 
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Figure  21.  Power  Spectral  Density  Curve  for  QPSK  Modem  With  Symbol  Rate,  600  SPS 
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Figure  22.  Power  Spectral  Density  Curve  for  QPSK  Modem  With  Symbol  Rate, 


Figure  23.  Power  Spectral  Density  Curve  for  BFSK  Modem  With  Data  Rate,  150  BPS 
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Figure  24.  Power  Spectral  Density  Curve  for  BFSK  Modem  With  Data  Rate,  300  BPS 
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Figure  26.  Power  Spectral  Density  Curve  for  BFSK  Modem  With  Data  Rate,  1200  BPS 
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Figure  28.  Power  Spectral  Density  Curve  for  BASK  Modem  With  Data  Rate,  150  BPS 
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BASK  MODEM 


Figure  30.  Power  Spectral  Density  Curve  for  BASK  Modem  With  Data  Rate,  600  BPS 


BASK  MODEM 
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Figure  31-  Power  Spectral  Density  Curve  for  BASK  Modem  With  Data  Rate,  1200  BPS 


Vita 


Robert  Leo  Jones  was  bom  August  19,  1947  in  Louisville, 
Kentucky.  He  graduated  from  high  school  in  Flaherty,  Ken- 
tucky in  June  1965  and  attended  the  University  of  Kentucky 
from  which  he  received  a degree  of  Bachelor  of  Science  in 
Electrical  Engineering  in  May  1970.  Upon  graduation  he  was 
commissioned  a Second  Lieutenant  in  the  United  States  Air 
Force  through  the  Reserve  Officers  Training  Corps  (ROTC) . 
Upon  completion  of  the  Communications-Electronics  Officers 
course  at  Keesler  AFB  MS  in  January  1971,  Captain  Jones  was 
assigned  as  Chief,  Communications-Electronics  Branch, 

2020th  Communications,  Squadron,  Shaw  AFB  South  Carolina. 
Captain  Jones  then  served  as  Site  Commander  of  a remote 
radio  relay  site  in  Eastern  Turkey  from  January  1974  to 
1975  whereupon  he  was  assigned  as  the  Chief,  Tactical  Con- 
trol and  Landing  Systems  Branch,  1868  Facility  Checking 
Squadron,  Rhein  Main  AB  Germany.  Captain  Jones  entered  the 
School  of  Engineering,  Air  Force  Institute  of  Technology 
in  June  1977 • 

Permanent  address:  305  Floyd 

Lynn  Haven  FL  32444 


93 


SECURITY  CLASSIFICATION  OF  This  PAGE  (Wien  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 


12.  GOVT  ACCESSION  NO. 


GE/EE/78-30y 


«.  TITLE  (and  Subdf/e) 


HEAD  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 


5.  TYPE  OF  REPORT  & PERIOD  COVERED 


Modem  Modeling  and  Characterization 


MS  Thesis 


7.  author*-*; 

Robert  L.  Jones 
Captain  USAF 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Air  Force  Institute  of  Technology 
(AFIT/EN) 

Wright-Patterson  AFB,  Ohio  45433 


11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

HQ  AFCS/OA 
Scott  AFB  IL 


4.  MONITORING  AGENCY  NAME  & AODRESS(i/  different  from  Controlling  Office) 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBER^ a) 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4 WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

December  1978 


13.  NUMBER  OF  PAGES 

102 


15.  SECURITY  CLASS,  (of  thla  report) 

Unclassified 


15a.  DECLASSI  FI  CATION/ DOWN  GRADING 
SCHEDULE 


16.  DISTRIBUTION  STATEMENT  (of  thia  Report) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  the  abatract  entered  In  Block  20,  if  different  from  Report) 


18.  SUPPLEMENTARY  NOTES 


Approyedfor  public  release;  I AW  AFR  190-17 

\ 


JOSEPH,  P.  H±PPS,  (Maj>r,  USAF 
Director  of  Information  /~/f-7 


19.  KEY  WORDS  (Continue  on  reverae  aide  if  neceaaary  and  identify  by  block  number) 


20.  ABSTRACT  (Continue  on  reverae  aide  It  neceaaary  and  identify  by  block  number) 

A method  for  modem  modeling  and  characterization  is  presented. 
The  properties  of  Fourier  transforms  and  the  properties  of 
correlation  and  power  spectral  density  for  periodic  and  ran- 
dom functions  are  reviewed.  The  modems  modeled  are  character- 
ized by  the  theoretical  calculation  of  the  power  spectral 
density  and  autocorrelation  function  of  the  output  signal 
from  the  modem  under  ideal  conditions.  The  input  signal — * *■  ^ 
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